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Abstract. In this PhD tutorial article, we present experiments with quantum 
degenerate mixtures of fermionic and bosonic atoms in 3-dimensional optical lattices. 
This heteronuclear quantum gas mixture offers a wide range of possibilities for quantum 
simulation, implementation of condensed matter Hamiltonians, quantum chemistry and 
ultimately dense and quantum degenerate dipolar molecular samples. 

We show how quantum degenerate mixtures of ^''K and '^^Rb are created in the 
experiment. We analyze stages of evaporative cooling and show how a dynamic mean 
field collapse occurs during the final stage of the evaporation [T] as a result of attractive 
interactions. The particle numbers observed in our experiment have only been limited 
by this mean field collapse, resulting in an excellent starting point for our experiments. 

We explore magnetic field induced Feshbach resonances and demonstrate tuning of 
interactions [2] between and ^^Rb by means of heteronuclear Feshbach resonances. 
We observe both stable attractively and repulsively interacting mixtures. We analyze 
the mean field energy of the condensate and find qualitative agreement with a simple 
model. By making the interaction strong and attractive, we induce a mean field collapse 
of the mixture. For strong and repulsive interactions, we observe phase separation of 
the mixture. 

When loaded into a 3D optical lattice, a whole zoo of novel quantum phases has 
been predicted for Fermi-Bose mixtures. We present the first realization of Fermi-Bose 
mixtures in 3D optical lattices as a novel quantum many body system [5]. We study 
the phase coherence of the bosonic cloud in the 3D optical lattice as a function of 
the amount of fermionic atoms simultaneously trapped in the lattice. We observe a 
loss of phase coherence at much lower lattice depth than for a pure bosonic cloud and 
discuss possible theoretical scenarios including adiabatic processes, mean field Fermi- 
Bose Hubbard scenarios, and disorder-enhanced localization scenarios. 

After considering this many-body limit of mixtures in lattices, we show how 
fermionic heteronuclear Feshbach molecules can be created in the optical lattice [1] as 
a crucial step towards all ground state dense dipolar molecular samples. We develop 
rf association as a novel molecule association technique, measure the binding energy, 
lifetime, and association efficiency of the molecules. We develop a simple theoretical 
single channel model of the molecules trapped in the lattice [5, which gives excellent 
quantitative agreement with the experimental data. 
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1. Introduction 

Nature distinguishes between two fundamental types of particles: fermions and bosons, 
depending on the spin of the particle. Particles with half-integer spin are called fermions 
and obey Fermi-Dirac statistics. As a result of Fermi-Dirac statistics, in a system of 
indistinguishable particles, at most one particle can occupy a given quantum state. 
Particles with integer spin are called bosons and obey Bose-Einstein statistics. Any 
single-particle eigenstate of a physical system can be occupied by an arbitrary number 
of bosons. 

While all of the basic building blocks of an atom are fermions, an atom as a whole 
has bosonic or fermionic character depending on its total angular momentum. For 
a gas of atoms confined in an external potential, the quantum statistical properties 
of the atoms become important at ultralow temperatures where the thermal deBroglie 
wavelength of the constituents is on the order of the interparticle separation. For bosonic 
atoms, this leads to the onset of Bose-Einstein condensation as observed for the first 
time in 1995 in a gas of ^^Rb at JILA [6j, in "^^Na at MIT [7] and for '^Li at Rice 
University [8] in the special case of attractive interactions. For fermionic atoms, the 
onset of degeneracy is less spectacular due to the absence of a phase transition. In an 
ultracold spin-polarized fermionic gas, the appearance of a macroscopic Fermi sea has 
first been demonstrated at JILA in 1999 [9j. 

While the pioneering work on quantum degenerate gases revealed important 
quantum phenomena such as interference, superfluidity and nonlinear atom optics, 
recent years have seen spectacular progress in the realization of novel strongly correlated 
systems with ultracold quantum matter. Strong correlations are observed either 
when the interactions between the constituents become very strong (e. g. at Feshbach 
resonances) or when strong confinement imposes stringent boundary conditions (e. g. in 
the periodic potential of an optical lattice). 

Atomic systems offer a high degree of control of both the external confinement 
and the interactions between the constituents. The latter has become possible with the 
advent of Feshbach resonances flOi \TT\ which allow control of s-wave and even higher 
order scattering between atoms by means of external fields. Feshbach resonances have 
been the key to a series of ground breaking experiments. For two-component fermionic 
gases, they have allowed the exploration of the BCS-BEC crossover and demonstrated 
that fermions and bosons are not as far from one another as it may seem: A Bose-Einstein 
condensate of diatomic molecules made from two fermionic atoms can be continuously 
transformed into a BCS state of atomic Cooper pairs [121 1131 [IH IISI [TSJ [T7] . 

The "confinement-induced" approach to strongly correlated phases was first 
proposed for a gas of repulsively interacting bosons in 1998 by D. Jaksch and 
coworkers [18]. In particular, it was demonstrated that bosonic atoms loaded into 
an optical lattice are an ideal model system for the simulation of the Bose Hubbard 
Hamiltonian (see [19]) known from condensed matter physics, and it was predicted 
that the phase transition from a superfiuid to a Mott-insulating state can be induced 
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by merely increasing the lattice depth. The theoretical prediction, together with the 
experimental demonstration of |20], highlighted the potential of ultracold atoms in 
optical lattices for simulations of quantum many-body systems and for the realization 
of strongly-correlated systems. Optical lattices have since been the key to the 
observation of intriguing phenomena such as a Tonks-Girardeau gas of atoms in a one- 
dimensional geometry [211 ES] or a Kosterlitz-Thouless transition in a ID lattice (2D 
geometry) |23]. Two-body bound states (molecules) in homonuclear systems have been 
engineered pll EH] and evidence for fermionic superfluidity has been reported in a cloud 
of ^Li loaded into a 3D optical lattice [26] . 

As a completely new area in the field of ultracold quantum gases, multicomponent 
quantum gases in 3D optical lattices have recently attracted a lot of attention. In the case 
of mixtures of fermionic and bosonic atoms, the different quantum-statistical behavior 
of the components gives rise to fundamentally novel quantum many-body phases. In 
the extreme case of pairing of fermions with one or more bosons, a whole zoo of new 
quantum phases of these "composite fermions" has been predicted [27]. Fermi-Bose 
mixtures in 3D optical lattices may exhibit fermionic pairing which is mediated by the 
presence of bosonic atoms in full analogy to solid state superconductivity, and there 
are interesting connections to high-Tc superconductivity [28], [291 130] . Even before such 
"atom pairs" form, Fermi-Bose correlations are predicted to become manifest in polaron- 
related physics of fermions dressed by a bosonic cloud [31] and quantum percolation [32] . 
These phenomena are connected to disorder induced localization scenarios. In reduced 
dimensionality, phenomena such as charge-density waves [HI [30] and supersolids [33] 
are predicted to occur. 

Any atomic Fermi-Bose mixture is necessarily also a special case of a heteronuclear 
system, i. e. a two-component system where both nuclei have a different decomposition. 
In some cases, such as ^Li - ''Li, the resulting mass difference is small, whereas in 
others, such as ^Li - ^^^Cs, it is very large. The mass difference, as we shall see, 
gives rise to interesting features already in the case of harmonic trapping of a mixture, 
but it ultimately opens up an interesting and highly promising approach to dipolar 
gases, quantum computation and simulation: if the components of such a heteronuclear 
mixture are brought together to form a molecule in its internal (rovibrational and 
electronic) ground state, the mass difference gives rise to a permanent molecular dipole 
moment. The resulting dipolar interaction, together with the high density and ultracold 
temperature of the initial atomic samples, would pave the way for novel quantum 
gases with dipolar interactions and quantum computation with polar molecules [34] . 
In more exotic heteronuclear systems, these techniques could be used for precision 
measurements [35l [36l [37] . 

From the experimental point of view, there has been an impressive series of 
experiments on homonuclear systems which we have partly mentioned above, but 
experiments on heteronuclear systems in lattices, which Fermi-Bose mixtures are a 
special case of, have been scarce. By the end of 2005, the only experiment with Fermi- 
Bose mixtures in optical lattices has been reported by Ott and coworkers at LENS [38] . 
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In these experiments, the "insulating" behavior of a trapped ideal Fermi gas in a ID 
lattice has been compared to coUisionally induced transport of fermionic atoms in the 
presence of a bosonic cloud. 

In this tutorial article, we present experiments performed at the University of 
Hamburg as part of our PhD theses [391 SO]- The article is organized as follows: we 
start our discussion by giving a cartoon picture of harmonically trapped Fermi-Bose 
mixtures and review a simple Thomas-Fermi model for the density distributions. As a 
function of the heteronuclear interaction, we identify regimes of stable mixtures with 
attractive and repulsive interaction as well as regimes of collapse and phase separation 
for a large heteronuclear interaction strength [H]. 

We show how we create large harmonically trapped Fermi-Bose mixtures of 
and *^Rb in a magnetic trap. We analyze stages of evaporative cooling and identify 
signatures of a dynamical mean-field collapse of the mixture as a result of attractive 
interactions. We show that large particle numbers of 7- 10^ and 1.2 ■ 10^ ^''Rb atoms 
can be achieved [T], only limited by the aforementioned mean field collapse. 

We show how heteronuclear interactions can be tailored by means of Feshbach 
resonances. The tunability of interactions gives access to the full phase diagram of 
the harmonically trapped mixture and allows us to produce both stable attractive and 
repulsive mixtures and to observe phase separation and collapse [2]. 

As a novel quantum many-body system, we discuss properties of Fermi-Bose 
mixtures trapped in 3D optical lattices. We show how already a small admixture 
of fermionic atoms reduces the bosonic coherence at much lower lattice depths than 
for the pure bosonic (superfiuid - Mott insulator transition) case [3]. We discuss 
various theoretical scenarios, from mean field models over thermodynamic processes 
and disorder-enhanced scenarios. 

Combining the ability to load the mixture into a 3D lattice and tunability of 
interactions, we demonstrate creation of heteronuclear Feshbach molecules in a 3D 
optical lattice [HE]. When combined with coherent Raman de-excitation schemes for 
the molecular ro- vibrational manifold, this constitutes a key step towards the production 
of all-ground state samples of dense, polar molecules. 

2. Fermions and bosons in a harmonic trap 

We will start our discussion with a basic introduction to density distributions of trapped 
interacting Fermi-Bose mixtures. This will help us identify regimes of stable mixtures 
as well as interaction-induced instabilities and give us a basic understanding of the 
phase diagram. Harmonically trapped interacting Fermi-Bose mixtures have first been 
considered in a paper by M0lmer in 1998 j41j, see also ref. [42]. For the bosons, we 
will neglect the influence of the thermal cloud and consider a Bose- Einstein Condensate 
in the Thomas-Fermi approximation. For the fermions, we assume the Thomas-Fermi 
approximation at T = 0. In this case, the bosonic and fermionic density distributions 
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nsir) and np{f) are given by the self-consistent solution to the following set of equations: 

nB{r) = ^— max [/is - gpB ■ npir) - VB{r), 0] (1) 
9bb 

npir) = max [fip - Vpif) - Qfb ■ nB{f),Of^'^ 

The Bose and Bose-Fermi interaction strengths are related to the corresponding 
scattering lengths according to 

gBB = 2n}faBB/lJ'BB (2) 

gpB = 2-Kh^aFB/ ^-FB (3) 

where ^bb and iipB are the reduced masses associated with the two-boson and the 
boson-fermion case. The bosonic and fermionic chemical potentials /i^ and are fixed 
by the particle numbers according to: 

Nb= nBifiB,x)d^x (4) 



POO 

Np = / np{fip,x)d^x 
Jo 



VB{r) and Vp{r) is the external trapping potential acting on bosons and fermions, 
respectively. In this treatment, the two components influence each other's density 
distributions through a mean-field potential given by 

Up^B{r) = gpB ■ np{r) (5) 
UB^pir) = gpB ■ nB{r) (6) 

acting on bosons {Up^b) and fermions (Ub^f), respectively. Let us assume a harmonic 
trap with trapping frequencies cu^y ^ for the boson and uj^^y^^ for the fermion. For 
simplicity, let us assume that the trapping potential is harmonic and the same for 
both bosons and fermions. This situation is given for example for magnetic trapping of 
stretched states. For a far enough detuned dipole trap, this would also be valid. In this 
situation, the trap frequencies are related according to uj^ /uj^ = ^JniB/mK- We can 
introduce rescaled coordinates 



which have the dimension of the square root of an energy: r'^ = Vb = Vp] the coupled 
Thomas-Fermi problem becomes: 

nB{r) = ^— max [fis - Qfb ■ np{f) - Ol (8) 
9bb 

npir) = max Up -r - gpB ■ nB{r),0\ (9) 

Note that the problem only depends on f = \/x\ + x1 + x|. In the new set of 
coordinates, equations Q become 

Nb= / 47r2f2 . UBif) df (10) 



\mp J [Up 
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(a) Non-interacting mixture (b) weak attractive interaction {gpB / 9bb — ~4) 




(c) weak repulsive interaction {gpE/gBB = +4) (d) phase separation {gpE/gBB = +8 ) 



Figure 1. Density distributions for trapped Fernii-Bose mixtures calculated in the 
Thomas-Fermi model. Parameters are Nb = Np = 10^, and ljb = 27r • 50 Hz. Also, 
a-BB — 98.98ao 43j, corresponding to ^^Rb in |F = 2,mp — 2). 



/ 2 \ 3/2 2 

Nf= / 47r=^r2 ■ Upir) dr , (11) 



where we have used the spherical symmetry of the problem in the new coordinates. 
As one can see, the above problem can be formulated solely in terms of the potential 
energy of the mixture, and the trapping potential only enters the calculation through 
the geometric mean trapping frequency up = '^/uJiuJ2^- As we shall see, this is 
particularly relevant when comparing mixtures produced in different experiments and 
different harmonic traps. 

In the following, we are going to analyze solutions to this coupled mean field scenario 
for various heteronuclear interaction strength combinations. Typically, equations ^ 
and ^ are solved in an iterative scheme by initializing both bosonic and fermionic 
distributions with corresponding non-interacting distributions. In a first step, because 
the bosonic density is usually much higher than the fermionic density, equation ^ 
is used to update the fermionic density, taking into account the bosonic mean field 
potential. The fermionic chemical potential is determined according to equation (11). 
The scheme continues with updates to bosonic, fermionic, . . . density distributions until 



Heteronuclear quantum gas mixtures 



7 



sufficient convergence is obtained. 

Figure [l] shows solutions to equations ([T| for various scattering length combinations. 
We will limit the discussion to repulsive interactions between bosons {asB > 0) ffist. 



Let us start the discussion by looking at the non-interacting case depicted in fig. l{a] 



In this case, we have the fairly broad distribution of the fermionic cloud which is due to 
the Pauli exclusion principle preventing two atoms to come too close to each other. The 
BEC, on the other hand, shows the characteristic inverted parabolic distribution. As we 
turn on a weak attractive heteronuclear interaction (see fig. l(b)[ ), we observe that the 
shape of the BEC is almost unchanged, whereas the Fermi sea shows a peak in the center 
in the overlap region with the BEC. For relatively weak interaction, we can almost stop 
the iterative solution process to equations ([T]) after applying ^ once, considering only 
the effect of the BEC on the fermions. The high density of the BEC creates a mean field 
'dimple' in which the fermions are trapped. This leads to the characteristic increase in 
fermionic density in the overlap region. 

A similar picture arises for a relatively weak repulsive interaction between bosons 
and fermions (fig. l(c)[ ). In that case, the fermions get pushed out of the center of the 
harmonic trap by the BEC, and we observe a lowered fermionic density in the overlap 
region. 

As the Fermi-Bose interaction becomes stronger and repulsive, the fermions can 
become completely expelled from the overlap region with the condensate, and we observe 



the onset of phase separation as seen in fig. 1(d) For approximately equal particle 



numbers of bosons and fermions, phase separation occurs as a dense core of bosons 
surrounded by a dilute shell of fermions, although there can be more exotic scenarios. 

For very strong and attractive interaction, fig. [2]^a) shows subsequent steps in the 
iterative solution scheme to the coupled mean field problem. In this case, the presence 
of the condensate pulls a substantial amount of fermions into the center of the trap. 
In the limit of strong interactions, there is a sizeable back action onto the condensate, 
which in turn contracts even more, leading again to an increased fermionic distribution 
at the center of the trap. Beyond certain critical conditions, there will not be any self- 
consistent solution to equations ([T]). This is the regime of the mean field collapse of the 
mixture, in analogy to the mean field collapse observed for pure BECs with attractive 
interactions. If a stable mixture is brought into the regime of the mean field collapse 
by either increasing the number of particles or by increasing the interaction, the sample 
will collapse in a 3-body implosion, where the ever and ever increasing density at the 
center of the trap is accompanied by strong 3-body losses. We will later see how these 
losses can be observed in the experiment. 

To summarize the behavior of trapped mixtures, let us have a look at fig. |2](b), which 
shows a simple phase diagram of a trapped mixture. The phase diagram shows areas of 
stable mixtures for small values of the heteronuclear interaction; for large and attractive 
interaction, we observe the onset of collapse; for large and repulsive interaction, we 
observe the onset of phase separation. Note that the phase diagram is for an equilibrium 
situation at T = only; in particular, the presence of excitations can lead to collapse 
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Figure 2. (a) Iterations of the mean field model in the case of instability (collapse), 
(b) Phase diagram of a harmonically trapped Fermi-Bose mixture for Nb = 10^ and 
Nf = 5- 10^ 



earlier than expected from equilibrium mean field theory (see below). Also, the presence 
of a dense thermal cloud at T 7^ and the corresponding mean field confinement will 
lead to collapse earlier than predicted by our cartoon picture. 

So far, we have only talked about the case of repulsive interaction between bosons. 
In the case of attractive interaction [asB < 0), the condensate has a very strong tendency 
to collapse all by itself already for very small particle numbers. 



3. Harmonically trapped Fermi-Bose mixtures with large particle numbers 

After this discussion of basic properties of trapped Fermi-Bose mixtures, let us describe 
the basic experimental setup and route to producing trapped mixtures with large particle 
numbers. Figure [s] illustrates the basic concept of the experiment. A vapor of and 
®^Rb is produced within a rectangular glass cell (2D-M0T cell) at a pressure of between 
10"^ and 10~^ mbar using alkali metal dispensers [H]; in the case of ^°K, we use enriched 
sources 05] because of the low natural abundance of ^*^K. From this background vapor, 
a cold atomic beam is produced using two pairs of horizontal, orthogonal laser beams 
forming a two-dimensional magneto-optical trap |l6l HT] for both species. The beam 
is coUimated in the horizontal direction and directed along the vertical axis of the 
apparatus. The flux towards the lower part of the apparatus is enhanced by using 
a resonant pushing beam which is optimized to match the velocity distribution with 
the capture range of a 3D two-species magneto-optical trap located in the lower part 
of the apparatus in a second glass cell. This region is separated from the 2D-M0T 
region by a differential pumping stage, thereby achieving UHV conditions in the 3D- 
MOT region. From the cold atomic beam, atoms are accumulated in the 3D-M0T for 
typically ten seconds with final particle numbers on the order of 1 ■ 10^° for ^^Rb and 
2 • 10^ for (5 • 10^ in the presence of ^^Rb). In the case of ^°K, we use a dark SPOT 
scheme [Ml HH ED] to overcome losses due to light-assisted collisions between atoms. 
This scheme has been crucial for achieving large particle numbers of in the presence 
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Figure 3. Sketch of the basic experimental setup, excluding all laser systems, optics 
and electronics. In the 2D magneto-optical trap, a cold atomic beam is produced, 
sent through a differential pumping stage and recaptured in a 3D-M0T. Atoms are 
transferred into a magnetic trap, sympathetically cooled into the degenerate regime 
and then transferred into optical traps and optical lattices for further experiments. 



of '^^Rb in the combined MOT. 

After accumulation of atoms in the MOT and subsequent polarization gradient 
cooling, the atoms are prepared in a spin state appropriate for magnetic trapping. In 
our experiment, we use a mixture of in the \F = 9/2, mp = 9/2) state and ®^Rb 
in the \F = 2,mp = 2) state. Atoms in these states are prepared using light on 
the corresponding cycling transitions starting from ensembles in the F = 9/2 and 
F = 2 manifolds, respectively. We then switch on the decompressed magnetic trap 
with trapping frequencies of = = uf^ = 27r ■ 11 Hz to allow for almost 
mode-matched and efficient loading of the magnetic trap |5T]. The magnetic trap is 
then compressed adiabatically in the radial direction during 1.5 s to final frequencies of 
LJ^^ = LV^^ = 27r • 257 Hz, while remains constant. We typically prepare 3 • 10^ ®^Rb 
atoms and 2 ■ 10'' atoms in the magnetic trap. 

Inside the magnetic trap, we perform forced evaporative cooling of ^''Rb. At 
a given magnetic field, the energy difference between neighboring Zeeman states is 
SFKb/d^K — 9/4 higher for ^^Rb than for ^"^K; the evaporative cooling therefore 
predominantly removes ^^Rb, while is cooled in the thermal bath of ^''Rb. Any spin- 
polarized fermionic gas alone would suffer from Pauli-blocking of s-wave collisions |52] 
at the temperatures in the magnetic trap, making evaporative cooling of such a system 
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very inefficient. 

The ^"K-^^Rb system is characterized by a relatively large and attractive 
background interaction [53j, although there was some initial controversy on the exact 
interaction strength. The large attractive interaction favors efficient evaporative cooling 
but, on the other hand, make the system susceptible to the mean field collapse of the 
mixture discussed above. 

Let us now look at various stages of the evaporative cooling process as the ^^Rb 
component condenses and the interaction between and ^''Rb affects both density 
distributions. When the ^''Rb atoms are not completely evaporated, various regimes 
of mixtures are accessible, ranging from dense thermal ^''Rb clouds of 10^ ^^Rb atoms 
right at the phase transition point interacting with a moderately degenerate Fermi gas 
of 2 ■ 10^ atoms to deeply degenerate mixtures with almost pure condensates. We 
achieve > 1 ■ 10^ atoms in the condensate coexisting with 7 ■ 10^ ^"^K atoms, limited by 
the onset of the collapse discussed below. 

Throughout this discussion, the "^"K images were recorded with a time of flight Te 
of 3 to 5 ms. With an axial trap frequency of uj^ = 2tt ■ 16.6 Hz, this means that 
ujf ■ Te is on the order of 0.3 to 0.5, and with the relatively large size of the Fermi gas in 
the trap, we can still deduce information on the axial spatial distribution in the magnetic 
trap from the time of flight image. In the experiment, this has allowed the extraction 
of information about the outer regions of the Fermi gas where it is noninteracting and 
about the inner regions where it overlaps and interacts with the dense Bose-Einstein 
condensate. Note that the expansion of the samples is also affected by the attractive 
interaction between '^''K and ^^Rb; however, in a cigar-shaped trap, the main effect is 
on the radial expansion where the confinement is strong [SU [55] . 

The ®^Rb atoms, on the other hand, are imaged after a total time of fiight of 
20 ms. All of the experiments described in this chapter are based on imaging using 
the two-species detection technique (see Appendix D) which has allowed us to extract 
information about both species on the same CCD chip in one run. 

During the final stage of evaporative cooling, various regimes of interacting mixtures 
are observed. These different regimes become most evident when analyzing the axial 
ID profile of the interacting clouds (see fig. [4]). Here, we have analyzed in more detail 
three different stages. In (a), the Fermi cloud exhibits a small peak in the center on top 
of a flat proflle. This peak is due to the Fermi-Bose attraction increasing the fermionic 
density in the center of the trap where the BEC is formed. Once critical conditions for 
the onset of the collapse have been met as in (b), the collapse localized in the overlap 
region removes the overlapping fraction of the Fermi gas. The hole which is left behind 
can be reflUed from the outer regions of the trap, and, as we shall see, this refllling 
process can even lead to a second collapse of the mixture. The peaked structure which 
is shown in (c) is stable over relatively long timescales of several 100 ms once the total 
^"K number has been reduced to undercritical values; for comparison, we have included 
the result of a Fermi-Dirac flt to the experimental data which clearly fails to account 
for the interactions with the Bose-Einstein condensate in the center. In contrast to the 
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axial coordinate [arb. iin.] 



Figure 4. Density profiles of interacting ''°K-®^Rb mixtures in various stages of the 
evaporative cooling process [1 . Copyright (2006) by the American Physical Society. 



situation depicted in (a), which results from overcritical initial conditions and is on the 
edge to collapse, the condition of subfigure (c) is a signature of Fermi-Bose interactions 
in a stable mixture and corresponds to fig. 1(b) in our introductory discussion. 



Observation of a mean field induced collapse 

In the discussion of the density profiles, we have claimed that the appearance of the hole 
in the center of the trap (see fig. |4](b)) is due to the mean field collapse of the mixture. 
The distinction from localized loss processes in the presence of high densities is most 
apparent when looking at the associated time scales. In order to study the time scale of 
atom loss, we will discuss a situation where an only slightly overcritical mixture has been 
prepared. During the first few milliseconds of the hold time, where the rf knife is held at a 
fixed frequency, the condensate will still grow out of the thermal cloud and reach critical 
conditions after some initial delay. The result of the measurement is shown in fig. [5] 
The figure shows the particle number contained in the central part of the absorption 
image which helps enhance the visibility of the effect. Between 10 and 20 ms after the 
start of the measurement, a sudden drop reduces the integrated particle number in this 
area to about two thirds of its original value. The atom number then remains constant 
for some time and at 100 ms again undergoes a second sudden drop, which is due to 
a refilling of the cloud in the center from the outer regions and a subsequent second 
collapse. Looking at fig. |5j one might argue that this refilling should become apparent 
in the center particle number going up again after the initial collapse. There is no such 
evidence from the data. One must however take into account that the central atom 
number from the time of flight image only approximately reflects the particle number in 
the overlap region. The occurrence of excitations and strong unequilibrium conditions 
following the initial collapse may modify the expansion behavior and in particular lower 
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Figure 5. Sudden drop in atom number as a result of the initial collapse, followed by 
one revival of the collapse out of a non-equilibrium situation 1 . Inset taken from [56j . 
Copyright (2006) by the American Physical Society. 



the collapse conditions for the second occurrence of this mean-field instability. The 
message of fig. [5] is the step-like feature which shows a very sudden drop of the particle 
number characteristic for this rapidly contracting mean-field implosion. From fig. |5} we 
can constrain this timescale to below 10 ms, which is incompatible with pure 3-body 
loss in the absence of the collapse. The inset in fig. |5] shows results from numerical 
modelling of the mixture in the regime of mean field instability based on hydrodynamic 
equations [56]. The numerical simulations qualitatively predict the same behavior for 
the particle number as a function of time. Note that the simulations assume different 
trap parameters, resulting in the different time scales. Also, the sharpness of the steps 
is a function of the 3-body loss rate assumed in the simulation. 

One of the experimental challenges in measuring data as in fig. |5]is that over many 
runs, the initial conditions must be stable enough so that the time evolution can be 
well controlled. This is most important with respect to initial particle numbers and 
transfer from the MOT to the magnetic trap. When large overcritical mixtures are 
initially prepared, the steps shown in fig. [5] can no longer be resolved individually due to 
initial fluctuations and the very complex dynamics. Instead, in this limit the collapse is 
observed as a very rapid overall decay when measurements with independently prepared 
samples are visualized together in a way similar to fig. [5] It would be an interesting 
perspective to analyze the behavior of the cloud in the regime of instability using phase- 
contrast imaging which would allow several subsequent nondestructive measurements of 
the same sample. 

Stability analysis of Fermi-Bose mixtures 

During the above discussion, we have mentioned conditions of criticality of mixtures and 
critical particle number of the mean field collapse several times. Here, we summarize our 
experimental findings on stable and unstable mixtures and discuss the relation to the 
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value of the s-wave scattering length for collisions between K and Wo. Two criteria 
for observation of mean-field instabilities have been used here: 

• Particle number combinations where the overall decay of the mixture is much too 
fast for normal 3-body decay in connection with the appearance of the pronounced 
hole are considered as unstable. 

• The observation of the step-like drop as discussed in context with fig. [5] is a sign of 
instability. 

Other situations are identified as stable. The resulting atom number combinations are 
plotted in fig.[6j Stable atom number combinations are plotted in blue, and combinations 
found to be unstable in the experiment in red. Atom numbers observed as stable are e. 
g- 

Nf = 7-10^; Nb = 1.2- 10^; cj = 27r ■ 91 Hz (12) 

The error bars (c) result from an estimate of the ^''Rb and ^"^K atom number uncertainty 
of 20% and 30%, respectively. The uncertainty in the ^''K atom number is less important 
than for ^^Rb, since in the regime considered here the dependence of critical conditions on 
Np is much smaller than on Nb- Also shown in fig. [6] are particle number combinations 
from two other experiments, namely the one at LENS (a) and the JILA ^°K-^^Rb setup 
(b). The particle number combination (a) of 

A^i. = 2 ■ 10^ Nb = 1.5 ■ 10^; = 27r ■ 91 Hz (13) 

has been reported as critical for the onset of collapse in ref. [57], but is an order 
of magnitude smaller than the stable particle number combination observed in our 
experiment. Data point (b) is the stable particle number reported at JILA in ref. |58] 
without any sign of instability: 

Nf = 8- 10^ Nb = 1.8 ■ 10^; = 27r ■ 86 Hz (14) 

which raises the question why mixtures which are observed as stable in this experiment 
and by the JILA experiment have been reported as unstable at LENS. 

When comparing critical particle numbers of different experiments, a few words on 
the role of the trap parameters are appropriate. Conditions for instability are of course 
influenced by the trapping potential. In general, the stronger the external confinement, 
the lower the critical particle numbers or the critical interaction are. The scale for the 
onset of instability is set by the geometric mean of the trap frequencies ui which is 
given above for all three experiments. The mean harmonic trap frequency used in this 
experiment and at LENS are equal, whereas the JILA trap is just slightly shallower. 

Within the Thomas-Fermi approximation for bosons and fermions, we have shown 
in our introductory discussion that for fixed harmonic mean trap frequency, the aspect 
ratio does not influence the point of onset of collapse. This trap and the LENS trap 
are therefore directly comparable within this approximation, and the JILA trap is only 
slightly shallower and should therefore be slightly less subject to collapse than the two 
other traps for the same particle number. 
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Figure 6. Stability of "'"K-^^Rb mixtures [1^ as a function of bosonic and fermionic 
particle numbers Nb and Np with error bars (c). For comparison, we include stable 
(b) and unstable (a) particle number combinations observed in other experiments (see 
text). Copyright (2006) by the American Physical Society. 



The only difference between the LENS trap and the trap used here is the trap 
aspect ratio (our trap is shghtly more elongated). This may play a role in the following 
situations: 

• Breakdown of the Thomas-Fermi approximation. This may occur when 
going to very elongated quasi- ID geometries (which is not the case in any of the 
experiments discussed here). The breakdown of the Thomas- Fermi approximation 
in the radial direction would result from radial trap frequencies on the order of kHz. 
The general effect of such a geometry is a stabilization with respect to collapse. 
This stabilization is responsible e. g. for the stabilization of bright solitons in a 
narrow window before the onset of collapse [59j. Another reason may be that 
particle numbers become below, say, 10^. Neither of these scenarios applies to the 
experiments discussed here. 

• Gravitational sag. A trap which is radially more tight features a reduced 
differential gravitational sag between atoms of different mass compared to the 
isotropic situation. This differential sag is given by: 

-'{^-^) ^^'^ 
and vanishes for equal trap frequencies of both components. Since the differential 
gravitational sag reduces the spatial overlap of the components, some of the mean 
field energy is already "used up" in compensating for the gravitational sag. The 
larger gravitational sag in the Florence experiment thus results in the LENS trap 
becoming "less subject to collapse" for the same critical particle numbers, while the 
opposite is observed in the experiment. The role of the gravitational sag has also 
been discussed in [60] . 
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Figure 7. Aspect ratio of a condensate forming in a ^°K- 
the pure BEC expected aspect ratio of 1.3 (blue line) [SO] 
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Therefore, for the purpose of comparing the above-mentioned experimental results, the 
difference in trap aspect ratio does not play a role and cannot account for the order of 
magnitude difference in critical particle numbers. 

It has been pointed out by the LENS group |53 that with the very sensitive 
relationship between critical particle numbers and the Fermi-Bose interaction parameter 
apB, the observation of the collapse instability can be a very sensitive way to determine 
apB- This picture is true with respect to equilibrium conditions; however, the 
observation of instabilities during evaporative cooling, where excitations can never 
be completely avoided, complicates the picture, as we shall see. In this sense, the 
observation of a given particle number combination as stable imposes an upper limit on 
the interaction parameter. 

Excitations 

In our discussion of the implications of observed stability limits, we have mentioned 
several times that the presence of excitations in the sample during evaporative cooling 
means that an observed stability limit is always an upper limit on the interaction 
parameter. To illustrate this, fig. [7] shows the aspect ratio of the condensate before 
the onset of collapse in some of our experiments. The strong deviations (factor of ~ 10) 
from the equilibrium aspect ratio of 1.3 may be due to fast evaporation compared to 
the axial trap frequency, thereby creating strong axial excitations. 

These excitations can in turn locally increase the density and lead to critical 
conditions even for particle number combinations which are stable in an equilibrium 
situation. For example, a twofold increase of the Thomas-Fermi radius of the condensate 
due to excitations would lead to a fivefold increase in the central density of the 
condensate, leading to the same densities required for a an equilibrium collapse at 
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apB = — 200ao. The same is true for any experiment observing the collapse during 
the evaporation ramp: In the absence of more sophisticated theory taking into account 
the excitations, observation of certain particle number combinations as stable provides 
an upper limit on scattering parameters. 

When putting the observed experimental stability limit from this experiment (which 
produced the largest degenerate mixtures so far) into mean field theory, we obtain 
—apB < 284 oq as an upper limit on the interaction parameter. This upper limit 
coincides with the value deduced from the first ^°K-'^^Rb spectroscopy experiment at 
JILA and the corresponding resonance assignment. The identification of more Feshbach 
resonances and an upgraded collisional model performed at LENS as well as an improved 
magnetic field calibration has resulted in a scattering length of —215(10) ao [61], and 
the upper limit reported here is compatible with that value. Since the same argument 
is valid for any experiment, there is no contradiction between the observed stable 
particle number combinations and apB = —215(10) Oq. Rather, the above comparison 
demonstrates limits of a comparison between equilibrium mean field theory and the 
experimental observation during evaporation, where strong dynamics and excitations 
can occur. 

It is an important result from this experiment that large stable '*°K-*^Rb particle 
numbers can be produced, which was not clear a priori, and to have identified 
experimental signatures of the mean field collapse in Fermi-Bose mixtures, which have 
later been seen in experiments with tunable interactions, as we shall see now. 

4. Tuning of heteronuclear interactions 

So far, we have investigated the behavior of the mixture in the harmonic trap at a fixed 
background scattering length. In the early days of ultracold quantum g lot of 

interest in a particular atomic species was motivated by its static scattering properties; in 
the case of '^''Rb e. g., the moderately large background scattering length, combined with 
favorable collisional losses (and technical considerations such as cooling laser sources), is 
the main reason for the fact that the overwhelming majority of BEC experiments relies 
on ^''Rb these days. This picture changed dramatically with the advent of Feshbach 
resonances. 

A Feshbach resonance is a phenomenon originally introduced in the context of 
nuclear physics |62]. In cold atomic gases, these resonance phenomena have been 
discussed since 1992, especially for Cs |63l [M] which at the time seemed to be the most 
promising candidate for the observation of BEC in a dilute atomic gas. A Feshbach 
resonance occurs when the energy of the colliding atoms in an open channel resonantly 
couples to a quasibound molecular state in a closed channel. A Feshbach resonance 
allows tuning of the interaction properties of an atomic gas from large and repulsive to 
large and attractive attractive and vice versa. The enhanced interaction properties near 
the Feshbach resonances are often accompanied by enhanced inelastic collisions. These 
inelastic collisions provide an easy tool for the identification of resonances, but can at 
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Figure 8. Cartoon picture of a Feshbach resonance. Such resonances occur when the 
energy of two atoms scattering from a molecular potential ('open channel') is close 
to resonance with a quasibound state in a different molecular potential curve ('closed 
channel'). In general, both channels have different magnetic moments, allowing the 
bound state in the closed channel to be tuned relative to the open channel disscociation 
limit. This allows the s-wave scattering length (plotted here) between atoms in 
the open channel to be tuned from large and attractive to large and repulsive. In 
particular, when the bound molecular state is slightly below the dissociation threshold 
(immediately below the magnetic field reconance center position Bq, indicated by the 
vertical dashed line), s-wave interactions between atoms are repulsive (apB > 0), and 
weakly bound molecules can exist in the closed channel. When the molecular bound 
state is slightly above the resonance (right of the dashed line), no stable Feshbach 
molecules are formed, and the interaction between atoms is attractive (apB < 0). 

the same time affect the timescales for experiments in the vicinity of the resonances. 
Feshbach resonances have been suggested to be inducible via static electric fields ^5], rf 
radiation [66], via far-resonant or near-resonant laser fields [671 EH] or via homogeneous 
magnetic fields [63|, [Ml EH [70] . Nowadays, magnetically induced Feshbach resonances 
are a very versatile and common tool in cold atom laboratories. They occur when the 
magnetic moment of the bound state in the closed channel is different from the magnetic 
moment of the open channel and the energy of the bound state can thus be continuously 
tuned through the open channel dissociation limit (see fig. [s]). 

The first observation of Feshbach resonances in cold atom scattering has been 
reported almost simultaneously by P. Courteille and coworkers in *^Rb [10] and by S. 
Inouye and coworkers in a BEC of ^^Na [TT]. Tuning of interactions was demonstrated by 
Inouye and coworkers in the vicinity of the resonance by a static homogeneous magnetic 
field. Magnetically induced Feshbach resonances have later been reported in almost 
any alkah atom: ^Li [71] and ^Li [721 [73], ^'^Na tUj, [71], ^^Rb [TOl [75], «^Rb [76] 
and ^^^Cs [77]. Recently, optically induced Feshbach resonances have been reported 
in [781 [79]. 

The availability of an additional knob to tune the scattering length has had 
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important consequences for the development of the field. On the one hand, the fact that 
interaction parameters no longer have to be taken for granted has enabled new species 
to be condensed which would otherwise have been very difficult to nearly impossible, 
such as *^Rb [80] or ^^^Cs [ST]. Tuning of interactions is needed in order to achieve 
efficient rethermalization in this case. On the other hand, a large class of experiments 
exploring degenerate gases had been studying the properties of the system in terms of 
atom numbers or trap frequencies, whereas in many situations the parameter which is 
physically interesting is the interaction, and Feshbach resonances allow this knob to 
be turned just by changing the magnetic field. ^Li, for example, is a system with a 
strong attractive background interaction, allowing the production of BECs of only a few 
thousand atoms before the onset of collapse [8]. Phenomena like the collapse have first 
been observed by varying the atom number in the sample [82], thereby only indirectly 
affecting interactions through the density. Feshbach resonances allow these phenomena 
to be directly observed as a function of interactions [83] . 

Feshbach resonances have allowed to enter the regime of strongly correlated systems 
and have made a whole class of experiments possible which have no analogon in 
neighboring fields such as solid state physics. One prominent example is the BCS- 
BEC crossover in two-component Fermi gases interacting at a Feshbach resonance [121 
[131 [m [151 [16] . Apart from the BCS-BEC crossover, Feshbach resonances have been the 
key to success in ground breaking experiments such as solitons [THl [72], collapse [83] . 
molecule formation [82|85l[86l[87l|88llMll90l[9ll|92l[25]. 

In any case, the impressively large variety of experiments relying on Feshbach 
resonances have exclusively made use of homonuclear Feshbach resonances, that is 
Feshbach resonances arising between spin states of a single atomic species. As of early 
2006, the study of heteronuclear Feshbach resonances in collisions between certain spin 
states of two different atomic species (either bosonic or fermionic) was still in its infancy. 
Heteronuclear Feshbach resonances had been identified for the first time in Fermi-Bose 
mixtures of ^Li - ^^Na at MIT [93] and in ^"K-^^Rb at JILA [91]. The identification has 
been reported through increased inelastic collisions and atom loss at the resonance, but 
no tuning of interactions has been reported. 

The control of heteronuclear interactions offers fascinating perspectives in Fermi- 
Bose mixed systems and Fermi-Fermi systems of different mass: In Fermi-Fermi systems, 
heteronuclear Feshbach resonances would open the way to the exploration of the BCS- 
BEC crossover in a two-species fermion mixture with unequal mass [95l [96] • In Fermi- 
Bose systems, a different and complementary approach to fermionic superfluidity has 
been suggested in which the interaction between fermionic atoms is provided by bosonic 
atoms taking over the role of phonons in the solid state superconductor [971 EH l99l [29]. 
When loaded into three-dimensional optical lattices, a wealth of different phases with no 
equivalent in condensed matter systems have been predicted to exist [29l [271 HOOl llOlj . 
In addition, heteronuclear Feshbach molecules [H I102[ 1103] which have been produced 
for the flrst time within this work (see text) could be used as a starting point for the 
creation of polar molecules with novel anisotropic interactions. 
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Here we demonstrate the tunability of heteronuclear interactions in the vicinity of 
a Feshbach resonance arising between and ^^Rb [2] (see also [104]). To provide 
evidence for tuning of interactions and thus a variation of the scattering length around 
these heteronuclear resonances, the interaction energy of the Bose-Einstein condensate 
coexisting with a quantum degenerate Fermi gas is measured in the vicinity of an s-wave 
heteronuclear Feshbach resonance. Additional and complementary information is gained 
from a measurement of interaction effects of the Fermi gas on an expanding Bose-Einstein 
condensate and vice versa. As a first application of heteronuclear Feshbach resonances, 
the phase diagram for harmonically trapped Fermi-Bose mixtures ranging from collapse 
for large attractive interactions to phase separation in the case of strong repulsive 
interactions [H] has been studied. Tuning of interactions has been the prerequisite 
for the first formation of ultracold heteronuclear molecules as reported later in this 
article. 

Tuning of heteronuclear interactions in stable mixtures 

Evidence for tuning of interactions in the heteronuclear ^"K-^^Rb mixture is provided 
in two different and complementary ways: 

1. In a first experiment, we study tuning of interactions by a measurement of the mean- 
field energy of the Bose-Einstein condensate as a function of magnetic field in the 
vicinity of a Feshbach resonance. Due to the heteronuclear interaction, the Bose- 
Einstein condensate is confined in the combined potential of the external dipole 
trap and the heteronuclear mean- field potential (see our introductory discussion). 
The latter becomes stronger as the heteronuclear interaction increases. Hence, the 
effective confinement of the Bose-Einstein condensate which determines its mean- 
field energy is changed. A measurement of this interaction-dependent mean-field 
energy is performed by a sudden switch off of all confining potentials including the 
additional mean-field potential and an observation of the time-of-fiight expansion 
of the condensate. A sudden switch off of the mean-field potential is realized in 
a good approximation by suddenly switching off the Feshbach field, reducing the 
heteronuclear scattering length to its background value. A related study has been 
done in the very first demonstration of tuning of homonuclear interactions in a 
Bose-Einstein condensate of ^^Na [11] in the vicinity of a Feshbach resonance. 

2. In a second experiment, we study the influence of the heteronuclear interaction on 
the time-of-fiight expansion of the Bose-Einstein condensate and the Fermi gas. 
When the heteronuclear interaction is left on during time-of-fiight, the expansion 
of the two clouds is either slowed down due to attractive interaction or influenced 
by repulsive interactions. The study is performed by a sudden switch off of the 
external dipole trapping potential while the Feshbach field is left on during time of 
flight. 

For our studies, we have used the broadest of the available s-wave resonances which 
occurs between the |9/2, — 9/2) state and the ^'^Rb |1,1) state. The resonance 
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position is at i?o=546.8(l) G as determined by the transition between strong attractive 
and repulsive interactions (see below). At the resonance, the scattering length varies 
dispersively as a function of the external magnetic field |69] 



I AB . 

apB = flbg ■ ( 1 - ^ _^ I • (16) 



In equation (16), Bq is the resonance center, Obg = — 185ao jl04] is the background 
scattering length between bosons and fermions in the respective spin state and AB is 
the width of the resonance which in our case has been predicted to be -2.9 G |105j . 
Note that the width of the Feshbach resonance can be negative. In our case, this means 
that the scattering length will diverge to a — > +oo below the resonance {B < Bq) and 
a — s> — oo above the resonance {B > Bq). A sketch of the expected variation of the 
scattering length as a function of magnetic field for the ^°K-^^Rb s-wave resonance at 
B = 546.8(1) G is given in fig. g 

Experimental Procedure For the preparation of "Feshbach-resonant" mixtures, we start 
with a magnetically trapped thermal ensemble of ^''K and ^^Rb in the low-field seeking 
but non-resonant |9/2,9/2) and |2,2) states. Slightly before reaching the critical 
temperature for condensation of the bosonic ^^Rb ensemble, we ramp up the optical 
dipole potential of a crossed dipole trap within 500 ms. This novel "magic" trap 
(see Appendix B ) is designed to eliminate the differential gravitational sag between ^^Rb 
and by operating at the magic wavelength of 807.9 nm. In the combined optical and 
magnetic potential, a last step of rf-induced evaporation is performed within 1 s and a 
quantum degenerate mixture is prepared. While ramping down the magnetic trapping 
potential, care has to be taken not to induce spin flips in either or ^^Rb affecting 
the polarization of the atomic ensemble. To this end, an offset fleld of a few ten Gauss 
is applied throughout the ramp-down process. Residual heating of the cloud during the 
ramp down is subsequently compensated for by a last step of purely optical evaporation 
in the crossed dipole trap. We typically end up with a quantum degenerate mixture of 
5 ■ 10^ ^"K and up to 2 ■ 10^ ^''Rb atoms and no discernible bosonic thermal fraction. 

In a flnal step, the '^^K |9/2,9/2) and ^^Rb |2,2) mixture is converted into a 
|9/2,— 7/2) and |1,1) mixture using rf and microwave Landau-Zener transfer schemes 
(see Appendix A ). The flnal transfer to the resonant |9/2, —9/2) state in '^"K is performed 
at two flxed magnetic fleld values by an rf Landau-Zener transfer after the fleld has 
settled. For studies above the resonance [B > Bq), where the heteronuclear interaction 
is expected to be attractive (apB < 0), the flnal transfer is performed at a magnetic fleld 
of 550.5 G. For studies below the resonance center {{Bq — AB) < B < Bq), where the 
heteronuclear interaction is expected to be repulsive (ofb > 0), the transfer is done at 
B = 543.8 G. At these magnetic fleld values, the heteronuclear scattering length is still 
moderately attractive in both cases and not too different from the background scattering 
length between and ^''Rb in the |9/2, —7/2) and |1, 1) mixture. This is important in 
order not to induce oscillations in both clouds due to the sudden change of the scattering 
length. After the flnal transfer at these two magnetic fleld values, the fleld is ramped to 
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varying values in the immediate vicinity of the resonance either from above or below. 
The ramp is performed within 50 ms, thereby changing apB adiabatically with respect 
to any other time scale in the experiment: apB/o-FB ^ i-^i so that the clouds remain in 
equilibrium within the field ram pIU [m [M]. We then probe the expansion of the cloud 
at different values of the heteronuclear scattering length and under different conditions. 



Evidence for tuning I: Mean-field energy of the BEC While the kinetic energy of a 
Bose- Einstein condensate in the trap is negligible (Thomas- Fermi limit), the energy of 
the BEC is dominated by its interaction energy given by 

Ej/N = — aBB {nB) ■ (17) 
m 

In this equation, qbb denotes the Bose-Bose scattering length, {ub) is the average density 
of the condensate, N the number of atoms of mass m. This mean-field energy of the 
condensate can be measured by time-of- flight expansion. During time-of-flight, the 
stored interaction energy is converted into the kinetic energy of a freely expanding 
condensate. 

At first glance, the interaction energy of the condensate, proportional to the mean 
density and the bosonic scattering length, seems to be independent of the value of the 
heteronuclear scattering length. However, the effect of heteronuclear interaction enters 
through the density of the Bose-Einstein condensate which is affected by the presence 
of the fermionic cloud. The trapped fermions, along with the heteronuclear interaction, 
act as an additional heteronuclear mean-field potential Up^Bir) (equation (|5|) forming 
an effective potential composed of the dipole trapping potential Vb and the mean-field 
potential Uf^b^'t) and determines the density distribution of the bosonic cloud ^^(r) 
via the Thomas Fermi equation (see equation ([T])). 

In the case of attractive heteronuclear interactions {apB < 0), the density 
distribution of the fermionic cloud interacting with the BEC is increased compared 
to the non-interacting case. The effective confinement of the BEC is thus increased 
with increasing attractive heteronuclear interaction in the vicinity of the resonance, 
enhancing the density and the mean- field energy according to equation (17). Along with 
the intrinsic repulsion of the Bose-Einstein condensate, characterized by the positive 
bosonic scattering length cbb = 100.4(l)ao |43j, the effect is observable through the 
increasing kinetic energy of the bosonic cloud after time-of-fiight. 

In the case of repulsive interactions, the condensate acts as a repulsive bump in 
the center of the fermionic cloud, pushing the fermionic cloud out of the center of the 
trap. The fermionic density overlapping with the center of the BEC will be reduced, 
whereas the fermionic density in the border areas of the potential increases. In this 
case, the additional fermionic curvature acting on the condensate will push the bosons 

II In case of a thermal ensemble, one would also have to require apB/apB <C 7ei, where 7e; is the elastic 
collision rate. In our experiment, we obtain uji < jei and apB/tiFB ^ is the stronger requirement. 
Note that in the immediate vicinity of the resonance, the adiabaticity condition may break down as 
a/ a becomes comparable to Wi 
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towards the center of the trap, thereby again increasing density and mean field energy 
of the condensate. However, the effect is smaller than the enhancement of density in 
the attractive case. 



Some numerical results The above intuitive picture can be supported by numerical 
simulations. The simulations rely on the mean-field model summarized at the 
beginning of this article (see |1T]) and provide a self-consistent calculation of the 
densities of the harmonically trapped fermionic and bosonic cloud, considering the 
additional heteronuclear mean-field potentials Uf-,b and Ub-^f- Figure |9]^a) shows 
the theoretically expected variation of the mean-field energy Ei of a condensate of 
10^ bosons in the vicinity of the Feshbach resonance at 546.8(1) G, interacting with a 
fermionic cloud of 5 ■ 10^ atoms via a varying heteronuclear scattering length aFB(-B). 
In the figure, we have plotted the interaction energy of the condensate as a function of 
B — Bq. The mean field energy has been calculated according to equation (17) and is 
normalized to the mean field energy of a pure condensate, i. e. equation (17), with n-Q 
given by the Thomas-Fermi profile of a pure condensate. The behavior of the scattering 
length in the vicinity of the resonance is determined by equation (16) with an assumed 
width of the resonance of AB = —2.9 G. We have plotted the interaction energy only for 
stable interacting mixtures. In the same diagram, a sketch of the expected heteronuclear 
scattering length is plotted. 

On the repulsive side of the resonance (low-field side), the interaction energy of 
the condensate increases by 20% compared to the non-interacting case. At the same 
time, the heteronuclear scattering length varies between at ~ —3 G and ~ 400 Oq 
before the onset of phase separation. On the high field side of the resonance, where 
heteronuclear interactions are attractive, the effect of heteronuclear interactions is much 
more pronounced: the mean-field energy of the condensate increases by almost a factor 
of two while the scattering length is increased by approximately the same factor. The 
effect of heteronuclear interactions on the mean-field energy of the condensate is thus 
strongest on the attractive side of the resonance whereas only a slight effect is observed 
on the repulsive side. 



Experimental data In the above discussion, we have seen that the mean field energy 
of a Bose-Einstein condensate in a Fermi-Bose mixture does depend on the value of 
the heteronuclear scattering length, although the dependence can only be calculated 
numerically. The interaction energy can be measured in a time-of-fiight expansion 
experiment where all confining potentials are suddenly switched off, thus converting 
the stored mean-field energy into kinetic energy. The kinetic energy per particle is 
just given by Ek/N = ^mv'^j^^ where Vrms is the mean-root-squared velocity. Both the 
atom number and the velocity Vrms can be extracted from time-of fiight images, more 
precisely frms from the size of the cloud after time of fiight. 

Figure ^h) shows a measure of the condensates mean-field energy per particle as 
a function of the magnetic field in the vicinity of the 546.8(1) G heteronuclear Feshbach 
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Figure 9. (a) Theoretically predicted variation of the mean-field energy of the 
condensate as a function of B — Bq in the vicinity of the 546.8(1) G Feshbach resonance. 
The mean-field energy is normalized to the interaction energy of a pure bosonic 
condensate Ej. (b) Interaction energy of a BEG per particle in the vicinity of a 
heteronuclear resonance normalized to the interaction energy of a pure condensate. 
On both sides of the resonance, the interaction energy increases although the effect 
is much more pronounced on the attractive side of the resonance where both the 
fermionic and the bosonic density distribution is enhanced. On the repulsive side 
of the resonance, the interaction energy increases only slightly. In the same diagram, 
we have plotted the expected variation of the heteronuclear scattering length across 



the resonance according to equation (16 1 



resonance. The experimental procedure for the preparation of resonant mixtures in the 
vicinity of the resonance has been described above. After the final magnetic field ramp 
from either above or below the resonance center to final field values near the resonance 
and a 10 ms hold time, both the Feshbach field and the external optical potential are 
suddenly switched off and and *^Rb are detected after time-of-flight. 

On the low-field side of the resonance, the interaction energy increases only slightly 
with increasing repulsive heteronuclear interaction. Within the investigated magnetic 
field range below the resonance, the scattering length is expected to vary from to 
-1-1000 flo- In fhfi same magnetic field interval, the interaction energy is observed 
to increase only by 25%. A comparable effect has been predicted by the numerical 
simulations discussed above. On the high-field side of the resonance, the effect of 
heteronuclear interactions is much more pronounced. The heteronuclear scattering 
length is expected to vary between — SOOoq and — lOOOao within the investigated 
magnetic field range. As can be seen from the data, the mean-field energy per particle 
stored in the condensate increases in the same interval by a factor of 3. This is in 
qualitative agreement with theory which predicts a pronounced increase in interaction 
energy of the condensate by a factor of 2. The remaining quantitative disagreement 
between theory and experiment might be due to a finite magnetic field settling time 
after the adiabatic ramps. 

For magnetic field values above B > 551 G, the heteronuclear scattering length 
is expected to recover slowly to its background value of — 185ao. Consequently, the 
interaction energy per particle is expected to saturate somewhat above the interaction 
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energy of a pure condensate. On the low magnetic field side, however, the scattering 
length changes sign from repulsive interaction to attractive a.t B ^ 543.9 G. Below this 
zero-crossing, the scattering length will become negative and eventually recover to the 
negative background interaction. Hence, the interaction energy per particle is expected 
to increase again below the zero-crossing converging towards the same interaction energy 
as far above the resonance. 



Evidence for tuning II: Interaction effects during expansion In a second experiment, we 
have studied the influence of the heteronuclear interaction on the time-of-flight expansion 
of the Bose-Einstein condensate and the Fermi gas. When the heteronuclear interaction 
is left on during time-of-flight, the expansion of the two clouds is either slowed down due 
to attractive interactions or influenced by repulsive interactions. The external dipole 
potential is suddenly switched off, whereas the Feshbach field is left on during time-of- 
flight expansion. 

The effect of attractive interactions on the time-of-flight expansion of the bosonic 
and fermionic cloud has already been studied in [54] in a Fermi-Bose mixture of and 
®^Rb. However, in these experiments, the heteronuclear interaction was determined by 
the background interaction strength and could not be tuned. The observations in [51] 
are therefore restricted to the observation of an affected time-of-flight expansion of the 
bosonic and fermionic cloud, respectively, in comparison to the expansion of a pure 
bosonic or fermionic cloud. 



Figure 10 shows the width of the bosonic cloud (a) and fermionic cloud (b) after 
7.5 ms and 25.2 ms time-of-flight, respectively. To account for enhanced particle loss in 
the vicinity of the resonance, the observed width data should be normalized by N^^^ 
and N^^^ for the bosonic and fermionic cloud, respectivelyj^ However, the observed 
particle number loss is on the order of 20% in the studied magnetic field range and the 
correction due to particle number variations is therefore below 5%. On the repulsive 
side of the resonance, where the effect of interaction on the time-of-flight expansion is 
very weak compared to the attractive case, particle loss is even smaller. In the same 
diagram, we have plotted the expected variation of the heteronuclear scattering length 
across the resonance. 

On the high-field side of the resonance, where heteronuclear interactions are 
attractive, we observe the width of both the fermionic and the bosonic cloud to decrease 
when approaching the resonance. This is a clear signature for increasing heteronuclear 
attraction. The simultaneous expansion of both clouds is slowed down due to the 
attractive potential that is present throughout the expansion process. In principle, the 
effect of the heteronuclear attractive interaction on the expansion of the BEG and the 
Fermi gas is twofold. On the one hand, the densities of both components are significantly 
increased inside the trap. This enhancement of density which corresponds to a tighter 
confinement of either of the clouds, will lead to a faster expansion (as discussed for the 

^ The equilibrium width of a non-interacting Fermi gas grows with N^^^ whereas the width of a Bose- 
Einstein condensate increases with iV^/^ 
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bosonic cloud above). On the other hand, the condensate and degenerate Fermi gas will 
interact during expansion and most importantly in the early stages of time of flight, 
thereby reducing significantly the expansion rate. This slow-down is determined by the 
relative dynamics between the condensate and the degenerate Fermi gas and dominates 
the experiments. 

On the repulsive side of the resonance, however, the effects of heteronuclear 
interaction are less pronounced. The width of the condensate after time-of-flight stays 
almost constant, whereas the width of the Fermi gas after time of-flight increases slightly. 
In the immediate vicinity of the resonance, this increase exhibits a change of slope. When 
the fermionic density at the trap center vanishes at even higher repulsion, the potential 
felt by the Bose cloud will rather be that of the pure external trapping potential with 
quite a sharp transition to a very steep higher order potential created by the fermionic 
density in the outer shell, at the edges of the condensate. We identify this region with 
the regime from 546.4 G to the center of the resonance at 546.8 G where, as seen in 
fig. ^h), the width of the condensate saturates. This may indicate that at complete 
phase separation, the repulsive interaction leads to a rapid expansion of the Fermi gas 
suddenly accelerated outside when the external potential is switched off and the repulsive 
bump of the BEG in the center maintained. 



Observing collapse and phase separation 

As we have seen so far, interactions between fermions and bosons can fundamentally 
affect the properties of the system such as the density of the components and expansion 
properties. However, when increasing the scattering length above a certain critical 
interaction strength, two instabilities can occur: collapse for the case of attractive 
interactions and phase separation for the case of repulsive interactions. In our 
experiments, we have been able to access the complete phase diagram of the mixture and 
observe clear signatures of both collapse and phase separation. The phase diagram has 
been accessed through tuning of interactions in the vicinity of the Feshbach resonance. 
Although the mean-field collapse of a harmonically trapped Fermi-Bose mixture has been 
studied before, in those experiments critical conditions have been achieved by increasing 
particle numbers in the system, holding the value of the heteronuclear interaction fixed 
(see |107[ [H HQ] and the discussion above in this article), whereas in the experiments 
reported here, the collapse is induced by tuning of the heteronuclear scattering length. 



Collapse For attractive interactions above a certain critical value, the Fermi-Bose 
mixture is expected to become unstable with respect to collapse. In the vicinity of 
the resonance occurring at 546.8 G, attractive interactions are strongest immediately 
above the resonance center. In fig. 10, the corresponding region has been shaded in 
grey. This region immediately above the resonance center has been studied in detail in 



fig. 11 Figure |ll[a) shows the atom numbers of both the fermionic and the bosonic 
cloud as a function of magnetic field when adiabatically approaching the resonance from 
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Figure 10. Observed width of the components after time-of-flight expansion as a 
function of magnetic field when the resonant interaction is left on during expansion [5]. 
(a) Width of the bosonic cloud, (b) Width of the fermionic cloud, (c) Heteronuclear 
scattering length. The region shaded in grey indicates the magnetic field regions 
where instability with respect to collapse and phase separation have been observed 
(see below). The black vertical line marks the observed transition from attractive to 
repulsive interactions (resonance center). (1 pixel = 3.2 /im; time of fiight 25.2ms for 
^^Rb and 7.5 ms for ^''K). Copyright (2006) by the American Physical Society. 



above. In the vicinity of 547.4 G (at a detuning of about 0.6 G above the resonance), we 
observe a sudden drop in both the fermionic and the bosonic atom number. We loose 
almost 50% of the fermionic atoms and a significant fraction of the bosonic cloud. In 
contrast to previous work, reporting the onset of instability as a function of atom number 
(see preceding discussion and [1_07]), the collapse is now due to tuning of interactions 
above a certain critical interaction strength in an otherwise undercritical mixture. The 
precise magnetic field value at which the collapse happens is of course dependent on 
the initially prepared atom number. In the experiments, we observe the onset of the 
collapse at a magnetic field value corresponding to a theoretically predicted scattering 
length of -1000 ao. 

To gain some information about the time scale on which the collapse happens, we 
have repeated the experiment in a slightly varied version: In this case, we have stopped 
the magnetic field ramp at the observed critical magnetic field value and then varied the 
hold time of the mixture, thereby observing the collapse happen as a function of time. 
This is illustrated in fig. [ll|^b). We cannot fully exclude some residual variation of the 
scattering length due to settling of the magnetic field after the nominal stop of the ramp. 
The onset of the collapse, again visible as a sudden drop in atom number, is retarded by 
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Figure 11. Induced mean- field collapse of the mixture |2j. (a) Sudden drop of atom 
numbers for critical heteronuclear interactions, (b) The onset of collapse is retarded 
by a timescale given by the trap frequency, (c) Sample time of flight images showing 
the bimodal distribution in the fermionic component and the sudden collapse of the 
system in the overlap region with the BEC. (d) As the collapse happens, the sample 
is excited and heated, visible in the width of the fragments. After the collapse, the 
mixture finds a new equilibrium. Copyright (2006) by the American Physical Society. 



a timescale given roughly by the trap frequency and happens on a timescale <1 ms. The 
dynamics associated with the collapse can be seen in fig. [ll|(c) and (d). Figure [Tl](d) 
shows two time of fiight images of the fermionic cloud prior to the collapse and after 
the collapse has happened (when leaving the Feshbach field on during expansion) . Prior 
to the collapse, the density of the fermionic cloud is strongly enhanced. The fermionic 
density profile after time of fiight shows a bimodal distribution with a broad and dilute 
pedestal and a very dense center. The latter part of the fermionic cloud is trapped inside 
the mean-field potential of the BEC. When critical conditions for collapse are achieved, 
this central core suddenly disappears. During the collapse, the overlap region of the 
fermionic cloud with the BEC is destroyed by a three-body implosion, leaving behind 
only the broad dilute pedestal. The three-body implosion causes significant heating 
and excitation in the remaining sample, refiected in the evolution of the width of the 
fermionic cloud as the collapse happens — see fig. [TT|^d). After the collapse, the mixture 
assumes a new equilibrium through evaporation and particle loss which is refiected in 
the decreasing width of the fermionic cloud after the collapse has happened. 

Phase separation Phase separation due to repulsive interactions in a composite system 
of harmonically trapped fermions and bosons has been discussed in theory [HI 1108] , but 
never been explored in experiment. Tuning of heteronuclear interactions has enabled 
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Figure 12. Vertical position of the Fermi g function of magnetic field [2]. 

The grey shaded area indicates the assumed region of full phase separation, where the 
fermionic density is expected to vanish at the BEC core. Due to gravitational symmetry 
breaking, the fermions are pushed above the BEC, an effect amplified by time of flight 
expansion. The inset shows the corresponding fermionic density distribution where 
most of the density is concentrated in the upper part of the image. Copyright (2006) 
by the American Physical Society. 



US to enter the regime of repulsive heteronuclear interactions, where phase separation 
occurs. In the limit of vanishing differential shift due to gravity and for our experimental 
parameters, phase separation will occur as a shell of Fermions surrounding a dense BEC 
core (as we have seen in the introductory discussion). 

Inside a harmonic trap and in the presence of gravity, atoms experience a 
gravitational sag given by —g/u'^. For systems with different masses, such as the ^°K- 
^''Rb system, this will in general lead to a differential gravitational sag between the 
components, as the trap frequencies may be different. Although our optical dipole trap 
has been designed to eliminate the differential gravitational sag, a slightly different 
gravitational sag has still been present in these experiments and breaks the symmetry 
of the system which therefore favors phase separation to occur in the vertical direction. 
As a consequence, the position of the fermionic component in the time of flight image 
is shifted upwards as a function of detuning from resonance, with an even stronger 
slope in the region of complete phase separation (see fig. 12). An important aspect is 
that the shift in position between fermions and bosons in the trap is amplified by the 
repulsive interaction during expansion if we leave on the interaction. Thus, the small 
initial symmetry-breaking in the direction of gravity is strongly enhanced and clearly 
visible in absorption images such as in the inset of fig. 12 , where the Fermionic density 
is concentrated in the upper part of the image. 
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5. Fermi-Bose mixtures in optical lattices 

Quantum degenerate gases in optical lattices provide unique opportunities for 
engineering of many-body Hamiltonians discussed in the context of condensed matter 
systems and even beyond in a very pure and controllable environment: The periodic 
potentials realized by optical lattices are almost defect-free. By tuning the depth of 
the lattice, the strength of the crystal can be controlled, thereby controlling e. g. 
tunnelling with respect to onsite interaction. Playing with lattices of different geometries 
and dimensionalities, one, two- and three-dimensional geometries can be realized and 
continuously transformed from one into another. In addition, the use of Feshbach 
resonances allows the precise control of interactions between the atoms even allowing 
for continuous tuning from repulsive interactions to attractive interactions and vice 
versa [63], [Ml EH [70], [TO] [H]. The possibihty of realizing strongly correlated systems 
with ultracold atoms loaded into 3D optical lattices has been pointed out in a seminal 
paper by D. Jaksch and coworkers in 1998 [18] who discussed the realization of the Bose 
Hubbard Hamiltonian [19] with ultracold bosons loaded into 3D optical lattices. The 
experimental observation of the quantum phase transition from a superfiuid to a Mott- 
Insulator state with repulsively interacting bosons followed in a pioneering experiment 
by M. Greiner and coworkers in 2002 [20]. Following experiments with either bosonic or 
fermionic atoms in optical lattices demonstrated the enormous degree of control inherent 
in these intriguing physical systems (see review |101] and references therein). 

In the introduction, we have mentioned the bright possibilities for heteronuclear 
quantum gases in 3D optical lattices as a novel quantum many-body system. In a 
heteronuclear Fermi-Bose mixture, the different quantum statistical behavior of the 
two components are predicted to give rise to fundamentally novel quantum many-body 
phases (see introduction). Here, we report on the first realization of a Fermi-Bose 
quantum many-body system confined in a 3D optical lattice [3]. In these studies, we 
investigate the coherence properties of bosonic atoms when interacting with a varying 
fraction of fermionic impurities. The interaction between fermions and bosons is held 
fixed at the attractive background interaction strength. When ramping up the optical 
lattice, we observe the phase coherence properties of the bosonic cloud to be strongly 
affected by a very small fraction [Nf/Nb ~ 3%) of fermionic atoms and find a fermion 
concentration dependent shift of the coherence properties with respect to the properties 
of a pure bosonic ensemble. A very small admixture of fermionic impurities induces a 
significant loss of coherence at much lower lattice depths as compared to the pure bosonic 
case. While the coherence properties of the pure bosonic system (loss of coherence with 
increasing lattice depth) can be explained in terms of the superfiuid to Mott-insulator 
transition and the associated many-body wavefunctions in the superfiuid and the Mott- 
insulating regime, the nature of the observed shift of the coherence properties towards 
lower optical lattice depth in the mixture is still a point of intense discussions in the 
community. Possible scenarios include thermodynamic effects like adiabatic heating 
when ramping up the optical lattice, disorder-induced localization scenarios or a shift of 
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the quantum critical point of the bosonic superfluid to Mott-insulator transition due to 
attractive interactions with the fermionic atoms (apB = —215(10) Oq [61]). In parallel 
to this work, similar studies and results have been reported at ETH Zurich |109j . 

5.1. Influence of fermions on bosonic coherence 

The coherence properties of a bosonic cloud can be revealed by a study of the bosonic 
interference pattern after a sudden switch off of the lattice and time-of-flight expansion. 
Whereas the observation of sharp interference peaks is a striking signature of long-range 
phase coherence of the bosonic cloud, decreasing interference peaks and the appearance 
of an incoherent background indicates the reduction of the coherence length. In a 
pure bosonic component, the loss of phase coherence with increasing lattice depth 
accompanies the transition from a superfluid to a Mott-insulating state. 

These experiments have been performed using a 3D optical lattice setup based on 
an Yb:YAG disc laser operating at a wavelength of 1030 nm. The lattice setup consists of 
three orthogonal, retro-reflected laser beams with orthogonal polarizations and a finite 
detuning of several MHz with respect to each other to minimize crosstalk between the 
different axes. Details of the lattice setup, including alignment, focussing, frequency 



and intensity stabilization can be found in Appendix C and |110[ IllH HOl 139] . 



In a first experiment, we have qualitatively studied the influence of a small 
admixture of fermionic atoms on the bosonic cloud. In these experiments, we have 
prepared 2 ■ 10^ fermionic atoms, coexisting with a pure Bose-Einstein condensate 
of 10^ ^^Rb atoms with no discernible thermal cloud in a combined magnetic and optical 
potential with trapping frequencies of u; = 27r ■ (50, 150, 150) Hz. Finally, the lattice is 
ramped up to various lattice depths between 2.5 Ef^ and 25 Ef-^. 

Figure [13] shows three sample images of the interference pattern of the pure bosonic 
cloud (top row) in comparison to the interference pattern of the bosonic cloud interacting 
with 2 ■ 10^ fermionic impurities (bottom row). The images have been taken after 15 ms 
time-of-fiight. To ensure that experimental conditions are really comparable in the 
two cases, we have prepared the pure bosonic component with the same experimental 
sequence as the mixture prior to ramping up the optical lattice: In both cases, we have 
prepared a deeply degenerate mixture of fermionic and bosonic atoms in the combined 
magnetic and optical potential. The pure bosonic cloud is then realized immediately 
before ramping up the optical lattice by a removal of the fermionic atoms from the trap 
with a short resonant light pulse, leaving behind a pure Bose Einstein condensate. The 



three sample images of fig. have been taken at lattice depths of 12.5 E^^, 20 Ef^ 
and 25 Ef'^, respectively. The lattice depths are given in units of the recoil energy for 
the ^^Rb component Ef-^ = {h^k'^) / {2m^^,) ^ h ■ 2.14 kHz where k = 2tt/X. 

In both cases, starting with a pure condensate or with a mixture of bosonic and 
fermionic atoms, we clearly observe a loss of interference contrast with increasing lattice 
depth marking the breakdown of long range phase coherence. As already mentioned, in 
case of the pure bosonic gas, the loss of coherence accompanies the well-known superfluid 
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Figure 13. Time of flight absorption images of the bosonic component 15 ms 
after switching off the lattice and trap potentials |3] (The colors and pseudo-3d 
representation encode the particle column density integrated along the imaging 
direction), (a) Pure bosonic ensemble of 10^ atoms at three different lattice depths, 
(b) Bosonic ensemble interacting with a fermionic impurity fraction of 20% for the 
same lattice depths as in (a) . We observe a striking loss of interference contrast in the 
images of the bosonic cloud interacting with fermions compared to the pure bosonic 
cloud. Copyright (2006) by the American Physical Society. 



to Mott-insulator phase transition [181 120l 1112] which occurs as a result of competition 
between the minimization of kinetic energy, parameterized by the tunneUing matrix 
element J which tends to delocalize the atomic wavefunction over the crystal and the 
minimization of interaction energy U (fig. [Tsj^a)). As can be clearly seen from fig. [T3|^b), 
the presence of fermionic impurities induces a loss of coherence at much lower lattice 
depths than for a pure BEC. At 12.5 E^"^, we observe sharp interference peaks in the 
pure bosonic cloud and the mixture. At 20 -E,^*^, however, the interference pattern 
vanishes almost completely in case of the mixed system, whereas significant interference 
is still visible in the case of pure ^^Rb. 

5.2. Characterizing the phase coherence properties of the bosonic cloud 

For a quantitative analysis of the observed phenomena, we adopt two different 
approaches appropriate for the characterization of the coherence properties of the 
bosonic cloud. These approaches have both been developed in the context of studies 
on the superfluid to Mott-insulator transition in an ensemble of repulsively interacting 
bosons: 

1. First, we characterize the interference pattern using the visibility of the interference 
fringes |113[ 1114] . The definition of visibility, in this context, is motivated by the 
definition of the interference contrast in optics. The visibility has been related to the 
correlation properties of the many-body wavefunction of the atomic cloud |1 13(11 14] . 
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(a) (b) 

Figure 14. Definition of the visibility of the bosonic interference pattern V |1131I114] . 
The visibility is calculated from the atom numbers in the first order interference fringes 
in comparison to the atom numbers in equivalent areas at intermediate positions 
between these maxima. The definition is illustrated by a typical interference pattern for 
shallow optical lattices (a) and with a sketch (b). Apart from the interference fringes, 
the sketch illustrates the first (light blue) and second Brillouin zone (dark blue). 



2. Second, we analyze the width of the central p = peak which is directly related to 
the inverse coherence length of the ensemble apart from some deviations due to the 
repulsive interactions between the bosonic atoms |115] . 

Visibility Figure [Mf^a) shows a typical interference pattern of a bosonic ensemble after 
time of flight released from a shallow optical lattice where a macroscopic fraction of the 
system occupies the lowest Bloch state and phase coherence is present across the lattice. 
An intuitive definition of the contrast of the interference pattern is given by |113[ 1114] 

^max '^min (18) 

^max ~l~ '^min 

where rimax is the total atom number in the first order interference fringes reflecting the 
p = 2hk momentum component and rimin is the sum of number of atoms in equivalent 
areas at intermediate positions between the maxima. 

At first sight, the relation of the visibility to the bosonic many-body wavefunction 
is not clear. However, it can be derived from the density distribution of the atoms after 
time of flight which is given by CIIl 

n{r) oc 

where w is the Fourier transformation of the Wannier function w{r) and 

S{k) = J2^Mik- in - r,)) (bib,) (20) 

is the quasimomentum distribution when we restrict ourselves to the lowest Bloch band. 
As S is the Fourier transformation of the correlation function (b\bj), sharp interference 
fringes are only possible when the correlation function varies slowly across the lattice. 
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Figure 15. (a) Visibility data for both a pure bosonic and a mixed system with 
an impurity fraction of 10%. We have fitted a phenomenological fit function to the 
data given by V(s) = ^"0/(1 + exp(Q! • (s — sq)). The initial visibility depends on the 
particular choice of the circles used to extract n,„ax and n„iin. (b) Width of the central 
interference fringe of the bosonic cloud both for a pure bosonic cloud and a mixed 
system with 10% impurity fraction. The width data shows a sudden increase above a 
certain lattice depth. The associated characteristic depth can be extracted from the 
data as the crossing point of two straight lines fitted to the data at low lattice depth 
and at high lattice depth. 



As outlined in |113j . this is the case of long-range phase coherence, \h\hjj extends only 
over a few lattice sites in the case of short range coherence, resulting in a reduction of 
contrast and visibility. Note that based on the definition of equation (18), the Wannier 
envelope |wp cancels out and the visibility V is directly related to the Fourier transform 
of the correlation function. 

Figure [l5|^a) shows the visibility of a bosonic cloud of 10^ atoms loaded into an 
optical lattice as a function of lattice depth. The data has been obtained starting from 
a pure condensate in the magic trap with trapping frequencies of u; = 27r ■ 50 Hz in the 
axial and radial direction. The optical lattice intensity is ramped up with a linear ramp 
of 0.5 £'^^/ms to various final lattice depths, nmax and nmin as defined in fig. 14 are 



extracted from the time-of-fiight images using a circle with a radius of 4 pixels around 
the center of the first order interference peaks and the displaced intermediate region, 
respectively. In the same figure, we have plotted the visibility of the bosonic cloud in a 
mixed system with 10% impurity fraction. 

Let us first study the visibility of the pure bosonic cloud. Starting from a lattice 
depth of 10 Ef°, we observe a continuous decrease in the visibility of the bosonic 
interference pattern which becomes particularly steep above 15 E^"^. Note that in our 
system, the quantum critical points for the onset of the superfluid to Mott-insulator 
transition in a pure bosonic system are lattice depths of 13.5 Ef-^ , 15.5 -E^"^, 16.9 Ef'", 
18.0 Ef^ and 18.8 Ef^ for n = l,...,n = 5, respectively. The observed loss of 
interference contrast is ascribed to the increasing fraction of atoms which enter the 
Mott-insulating phase. The quasi-continuous decrease in visibility, instead of a sharp 
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one above just one lattice depth corresponding to the quantum critical point of the 
transition, can be explained by the inhomogeneity of the system where atoms in different 
shells enter the Mott-insulating state with n atoms per site at different lattice depths. 



Phenomenological fit function For the analysis of the visibility data, we use a 
phenomenological fit function defined by 

= TT A ^ (21) 

1 + exp [a ■ [s - So)) 

where Vq is the initial visibility in the limit s ^ 0, Sq is a measure for the "critical" 
lattice depth and a is an additional fit parameter accounting for the steepness of the 
loss of contrast]^ As can be seen in fig. [T5|^a), the phenomenological function describes 
the experimental data quite well. In particular, the fit function accounts for the limiting 
cases of the superfluid state with a theoretical visibihty of V = 1 and the pure Mott-state 
with a visibility of V = 0. 



Comparison of pure bosonic to mixed systems Based on the phenomenological fit 



function of equation (21), we can now compare the decrease in visibility for a pure 
bosonic cloud to the decrease in case of a mixed system. The corresponding data is 
plotted in fig. [l5|^a). Each data point corresponds to 5 to 10 averaged measurements. 
The errorbars denote the standard deviation of the 5-10 measurements. To both the 
visibility data of the bosonic and the visibility data of the mixed system with 10% 



impurity fraction, we have fitted the phenomenological fit function of equation (21 ). The 



obtained fit parameters are summarized in table [T] Whereas Vq and the fit parameter 



System 




a 


Vo 


pure Bose cloud 


16.8 ±0.2 


0.45 ± 0.03 


0.8 ±0.02 


10% impurities 


13.7 ±0.3 


0.5 ±0.05 


0.8 ±0.02 



Table 1. Comparison of the parameters of the phenomenological function for a fit to 
two data sets, one for a pure bosonic system and one for a system with 10% impurity 
fraction. The errors correspond to the fit errors only. Note that e.g. the lattice depth 
calibration is estimated to be accurate within 5%, resulting in an error bar on sq of 

^ 1 



a are equal in the two cases within the fit uncertainty, we obtain a significantly lower 
"critical" lattice depth in the case of the mixed system than in the pure bosonic case. 
We observe a shift in the characteristic parameter sq of As = Sq^'^^ — s™^°'^ ~ — 3-E^^ 
towards lower lattice depth when comparing the mixed to the pure system. We have 
checked that we do not observe a significant variation in s™'^*^'^ for the mixed system 
when fixing a and Vo to the values of the pure bosonic system. 

For typical fit results with a « 0.4 and sq ~ 15, we obtain exp(— asp) « 0.002 and therefore 
V(s ^ 0) « Vo 
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Figure 16. Visibility data both for a pure bosonic and a mixed system with an 
impurity fraction of 10%. This time, we have plotted the visibility against U/zJ where 
U/J = ^ exp(2Vs) with s = Va/E^^. 



To estimate the reliability of the analysis based on the phenomenological fit 
function, we have analyzed the visibility data a second time, this time based on an 
analytically known dependence of the visibility in the limit of deep optical lattices. 
Based on perturbation theory in the limit of very weak tunnelling, the ground state of 
the system can be approximated by 

I^)^I*)mi + ^ E^1^^ I^)mi (22) 

<i,j> 

where is the Mott-insulating state. The analytical expression for the "column" 

integrated visibility is then given by [114j 

V oc ^ (23) 

where z is the number of nearest neighbors (z = 6 in case of a simple cubic 3D lattice). 
The above approximation is valid in the limit of deep optical lattices where tunnelling 
is strongly suppressed. 

To compare our analysis of the visibility data based on the phenomenological fit 
function to an analysis based on the known dependence of the visibility in the limit of 
deep optical lattices, we have analyzed the visibility both for the pure bosonic cloud and 
the mixture in the limit of deep optical lattices. 

Figure [16] shows the visibility data discussed above, this time plotted as a function 
of {U / zJ) which is directly related to the lattice depth in the tight-binding limit via 

l^^expPV^). (24, 

To extract a shift from the data, we have fitted a function C (^) in the limit of deep 
optical lattices {U/zJ > 10) and a constant in the limit of shallow optical lattices. 
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From these fits, we can again extract critical lattice depths for both the pure bosonic 
cloud and the mixture. Here, the critical depth is defined as the intersection point 
between the two fits. We observe a critical {U/zJ)c ~ 9.9 in the case of the pure 
bosonic cloud which corresponds to a lattice depth of 15.6 E^"^ and a critical depth of 
{U / zJ) PS 3.3 corresponding to Vq ~ 11.6 E^"^ in the case of the mixture. The exponent 
V is comparable to each other in both cases. We observe v = —1.3(1) for the pure 
bosonic cloud and u = —1.2(2) in case of the mixture. From the analysis, we extract a 
shift of ~ —4 Ef'^ which is comparable to the shift extracted from the phenomenological 
fit function. However, the extracted shift depends critically on the chosen fit interval 
for the deep lattice limit and varies easily by ±1 Ef'^, whereas the analysis based on 
the phenomenological fit function is independent of the fit interval. 

Width of the central interference fringe The phase coherence properties of the bosonic 
cloud loaded into the optical lattice can also be characterized by the width of the central 
interference fringe w which is directly related to the correlation length ( of the system via 
( oc w^^. The analysis is again motivated by fundamental properties of the superfiuid 
to Mott-insulator transition of a pure bosonic system. Whereas, in the superfiuid phase, 
the correlation length diverges, ( becomes finite in the Mott-insulating state and is on 
the order of the inverse energy gap of the system ( oc |115j . 

Figure [l5|(b) shows the half-width of the central interference fringe as a function 
of lattice depth for the pure bosonic and the mixed system for the same set of data as 
in the visibility analysis. In both cases, the width stays constant for some time or even 
slightly decreases for shallow optical lattices. Above a certain lattice depth, the width 
increases rapidly. The characteristic lattice depth for the sudden increase is extracted 
from the data by fitting two linear fits to the descending and ascending branches of the 
data. The characteristic lattice depth is then defined by the intersection point of the 
two linear curves. 



System 




pure Bose cloud 


16.5 ± 1.6 


10% impurities 


12.5 ± 1.5 



Table 2. Comparison of the "critical depth" of a pure bosonic cloud to a mixed system 
with 10% impurity fraction as extracted from the width data. 

Table |2] summarizes the characteristic lattice depths for the two cases. The error 
bars are again based on the fit uncertainties. Whereas the critical lattice depth of the 
pure bosonic cloud is given by 16.5(1.6) Ef-^, in the mixed system, we observe a critical 
depth of 12.5(1.5) Ef^^. The critical depth is thus shifted by -4.0(2.0) E^^ towards 
shallower optical lattices. The result is comparable to the respective results extracted 
from the visibility data. 
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5.3. Particle number dependence 



Based on the techniques of section 5.2, we have studied the change of coherence 
properties in the bosonic cloud as a function of particle number of the fermionic 
impurities. To this end, we have prepared mixtures of ~ 10^ bosons with a variable 
fraction of admixed fermionic atoms in the optical lattice. In the experiments, we have 
prepared mixtures with 3%, 7%, 10% and 20% fermionic atoms corresponding to 3 ■ 10^, 
7 ■ 10^, 1 ■ 10^ and 2 ■ 10^ fermions, respectively, in the combined magnetic and optical 
potential with trapping frequencies of a; = 27r ■ (50, 150, 150) Hz (for the experiments 
reported in [3]) or in the "magic trap" with trapping frequencies of = 27r ■ 50Hz. 
When choosing the above fractions of fermionic impurities, we have taken care to realize 
mixtures where the extension of the fermionic atoms is mostly still within the extension 
of the bosonic cloud in the optical lattice (details in [39]). 



Figure 17 summarizes our findings on the particle number dependence. We have 
analyzed samples with 10^ bosons with varying impurity concentration between 3% and 
20%. In case of the visibility data, the error bars in the diagram correspond to the 
"sample" rate of the underlying data. A sample rate of 1 — 1.25 Ef'^ has been translated 
into an assumed maximum error of the same size; this error is larger than the error 
based on the fit errors. In case of the width data, we have calculated the error bars 
based on the underlying fit errors of the involved parameters. As can be seen from 



fig. [17} an increasing impurity concentration leads to a considerable shift of the bosonic 
coherence properties towards shallower optical lattices. The shift increases with impurity 
concentration and reaches —5 Ef-^ for an impurity concentration of 20%. Note that the 
different data points correspond to ensembles that have been prepared on different days 
and with two different trap geometries. Whereas most of the measurements have been 
done with an external harmonic confinement given by a combined magnetic and optical 
potential with an aspect ratio of 1 : 3, some of the data points have been recorded 
with an almost spherically symmetric confinement of = 27r • 50 Hz of the 'magic trap' 
discussed in our experiments on tuning of interactions. In all these cases, we observe 
qualitatively and quantitatively the same characteristic behavior. 



5.4. Studies on systematic effects 

To exclude trivial irreversible heating effects as a source of the observed systematic shift 
in the coherence properties of the bosonic cloud, we have studied the reversibility of the 
observed coherence loss. To this end, we have ramped up the optical lattice to final 
lattice depths between 15 Ef-^ and 20 Ef-^, where the coherence of the bosonic cloud is 
either completely lost or at least significantly distorted, and have subsequently ramped 
down the lattice with the same ramp speed to varying lattice depths between 5 Ef'^ and 
7.5 Ef'^. Apart from some signs of non-adiabaticity visible in a slight asymmetry of the 
obtained visibility data with respect to the reversal point at deep optical lattices, we 
have been able to restore the initial visibility of the bosonic cloud almost completely. 
The maximum observed irreversible loss of visibility is on the order of 10% for both the 
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0.05 0.1 0.15 0.2 

fermionic impurity concentration 



Figure 17. Observed shift of the "critical" lattice depth characterizing the loss of 
coherence in the bosonic cloud as a function of fermionic impurity concentration. 



Fermi-Bose and the pure bosonic case. In case that the overall ramp sequence has not 
been fully reversible, allowing for an additional equilibration time of 5 — 10 ms after 
ramp-down has resulted in a near complete recovery of visibility. 

Note that the recovery of the visibility is not necessarily an indicator of adiabaticity, 
whereas an incomplete recovery would be a sure indicator. Also, there could be technical 
noises which have a more drastic effect on the mixture than on the pure bosonic system. 

5.5. Discussion 

There is currently an intense discussion on the nature of the observed shift of the 
coherence properties of the bosonic cloud due to a small admixture of fermionic atoms. 
Possible scenarios include a shift of the superfluid to Mott insulator transition of the 
bosonic cloud due to interactions with the fermionic impurities, thermodynamic effects 
such as adiabatic heating and effects of disorder induced localization scenarios. In the 
following section, we will briefly review the current status of the discussion and provide 
the interested reader with some simple estimates. 

Shift of the superfluid to Mott-insulator transition 

One of the possible scenarios is a shift of the superfluid to Mott insulator transition of the 
bosonic cloud due to the admixture of attractively interacting fermions. However, there 
is even a controversy about the sign of the expected shift depending on the underlying 
theoretical model: the Fermi-Bose Hubbard approximation or an effective potential 
approach. 

Fermi-Bose Hubbard approximation The change of the quantum critical point of 
the bosonic superfluid to Mott-insulator transition in a homogeneous system due to 
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the presence of fermionic atoms has been discussed in a recent work by Pollet and 
coworkers |119j . Based on the Fermi-Bose- Hubbard Hamiltonian |100j . Pollet et al. have 
discussed in zero and first order perturbation theory the change in the quantum critical 
point for the superfiuid to Mott insulator transition point {U /zJ)c due to interactions 
with the fermionic cloud. Here, we summarize the main results: Let us consider the 
Bose-Fermi Hubbard Hamiltonian of a homogeneous system 

Hbf = -JbJ2 ^l^^ + l^^^ E - 1) (25) 

<i,j> i 

-JfY. (26) 

<i,3> 
i 

(notation see |100j ). Under the influence of the fermions as a perturbation of the 
pure bosonic Hubbard Hamiltonian, the lowest order effect that the fermions have 
on the bosons can be derived when replacing the fermionic particle number operators 
nf by their expectation values. In this case, the perturbation to the pure bosonic 
Hubbard Hamiltonian is: UpB^^nf (^nf). To lowest order perturbation theory, the 
fermionic expectation values {nf^ are calculated by solving the pure fermionic Hubbard 
Hamiltonian (neglecting the effects of the bosons on the fermionic cloud). In the 
homogeneous system considered here, (nf) is independent of the lattice site (nf) = 
{n^) = {m). To lowest order perturbation theory, the fermions thus induce an overall 
shift of the system's energy and do not change the phase diagram of the bosonic cloud 
at all. 

The next order effect can be derived when adopting the results of linear response 
theory (see jll9 ]). In this case, one assumes the unknown bosonic density nB{q) to 
induce a fermionic density {np^q)) = UFB'^{T,q)nB{q) where E{T,q) is the Lindhard 
response function known from the theory of screening when calculating electronic charge 
distributions in condensed matter physics |120[ I12H I119j . The induced fermionic density 
distribution results in a back-action on the bosons and thus in an effective interaction 
for the pure bosonic cloud given by = Ubb + Up^EiT, q). As the Lindhard function 
is always negative, the interaction matrix element between two initially repulsively 
interacting bosons is effectively reduced by an amount Up^E(T,q) independent of the 
sign of UpB- This induced attractive interaction is in analogy to phonon induced 
attraction between electrons in conventional superconductors and is at the heart of 
recent proposals on boson-induced Cooper pairing of fermions in optical lattices |29j . 
Based on these considerations, the superfiuid to Mott-Insulator transition is shifted 
towards deeper optical lattices independent of the sign of Bose-Fermi interaction Ufb- 

The situation becomes more complicated when considering an inhomogeneous 
situation where the atoms experience both an optical lattice and an external harmonic 
trapping potential. In the case of attractive interactions, the equilibrium densities of 
both components will be enhanced in the center of the confining potential. However, 
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enhanced density results in an increased mean bosonic occupation number on a single 
lattice site, thereby shifting the transition from a superfluid to a Mott insulating state 
towards deeper optical lattices. 



Mean-field approach — Induced effective potential. Let us now come back to the 
scenario of a homogeneous optical lattice without additional overlapped harmonic 
potential. It has been suggested by Marcus Cramer |122] to consider the fermionic 
density distribution in the lattice as an additional mean-field potential U^p{r) which 
enters into the bosonic many-particle Hamiltonian via the confining external potential 
VeSectivc{r) = Vpcriodic (^) + U[jp{r). To lowcst ordcr perturbation theory, the additional 
fermionic mean-field potential is determined by the unperturbed fermionic density 
distribution given by np{r — r*j) = (m) \w^{r — rj)|, resulting in 

Umf = npi^UBF = (m) Ubf Y1 l^^i^' ^^\^ (27) 

i 

where is the fermionic Wannier function for the case of a pure sinusoidal external 
lattice potential of depth l^attice- Taking into account this additional potential, the 
bosonic many-particle Hamiltonian reads: 

Hb= J d^r¥j,{r} + V^periodic(rO + npirjUBF^ ^B(rO (28) 

+ ^Ubb J c/V^}3(f)^^(f)*B(r)^B(r) 

The fermions thus alter the properties of the pure external periodic potential and change 
both the depth and geometry of a single lattice well. From the effective many-body 



Hamiltonian of equation (28), we can derive an effective Bose- Hubbard Hamiltonian. 
The presence of the fermionic mean-field potential then changes both onsite interaction 
and tunnelling matrix elements of the bosonic atoms. 

In a recent numerical simulation, D. Liihmann and coworkers |123j have shown that 
the attractive interaction between and *^Rb of —215 oq will lead to a mutual trapping 
of particles. In particular, the bosonic effective potential will be strongly deformed and 
enhanced by the presence of the fermions and will thus result in an effective shift of the 
superfluid to Mott phase transition for the bosonic cloud. The authors derive a shift 
of the transition between 2.2 Ej. to 5.2 Ej. towards shallower lattices, depending on the 
filling factor of the bosons in the optical lattice and also observe a strong dependence 
on the interaction properties between bosons and fermions characterized by apB- Note 
that the main discussion in this paper assumes a homogeneous optical lattice. The case 
of an inhomogeneous optical lattice, i. e. an optical lattice with additional harmonic 
confinement in three dimensions (as realized in the experiment) is not covered. 



Adiabatic heating or cooling 

When adiabatically loading either bosonic or fermionic quantum gases into an optical 
lattice, the density of states for the components changes markedly, thereby e. g. changing 



Heteronuclear quantum gas mixtures 



41 



the characteristic temperatures for degeneracy in bosonic and fermionic quantum gases, 
Tc and Tp. Adiabatic loading of atoms occurs at constant entropy S, not at constant 
temperature, and the absolute temperature may change. The change in degeneracy 
characterized by T/Tc and T/Tp depends on the starting conditions such as temperature 
and atom number as well as the trap geometry and the interactions present in the system. 
Note that the temperature change is reversible as it occurs at constant entropy and is 
a purely thermodynamic effect. However, increasing T/Tc with increasing lattice depth 
may lead to a decreasing condensate fraction which might result in loss of coherence of 
the bosonic cloud. 

Adiabatic temperature change at constant entropy due to adiabatic loading of atoms 
into an optical lattice has been discussed in the literature for a homogeneous non- 
interacting Fermi gas |124] and for homogeneous and inhomogeneous interacting and 
non-interacting bosonic clouds |125[ 1126] . 

Recently, we have |127j analyzed thermodynamic properties of the adiabatic 
loading of an attractively interacting harmonically trapped Fermi-Bose system into an 
optical lattice. The analysis in the limit of vanishing optical lattice and deep optical 
lattice reveals that a significant adiabatic heating of the mixture is expected during 
the adiabatic loading process. While regions of adiabatic heating and cooling can 
be identified in M. Cramer's numerical simulations |127] depending on the starting 
conditions of the loading process, the presence of fermions leads either to a more distinct 
heating of the mixture or a less distinct cooling. In any case, the Fermi-Bose mixture 
will have a significantly higher temperature at the end of the loading process than the 
pure bosonic system. Although this temperature increase is a thermodynamic effect due 
to the adiabaticity of the loading process and can thus be reversed, it will definitely 
result in a decrease of the condensate fraction of the bosonic cloud in the optical lattice 
and thereby affect the phase coherence properties of the bosonic cloud in the lattice as 
characterized in the visibility and the width of the bosonic interference pattern. 



Disorder- enhanced localization scenarios 

Let us finally review intuitively possible disorder-related localization scenarios which 
might be observable in Fermi-Bose mixtures and would possibly also lead to reduced 
coherence in the bosonic cloud. For illustration purposes, let us start with a pure bosonic 
superfluid (fig. [Tsj^a)). Adding fermionic impurities and considering the attractive 
^°K-*^Rb Fermi-Bose interaction energy as an additional potential for the bosons, the 
"defects" caused by the fermionic impurities can be described by a local change of 
the effective optical lattice depth for the bosons due to the interparticle interaction 



(fig. 18 b)). If the energy level shift caused by the interaction energy is large enough, 
the superfluid bosonic wavefunction will not extend into this defect region, but will be 
scattered by the impurity. If scattering becomes frequent, interference effects along a 
closed scattering path are predicted to suppress transport and lead to a localization 
scenario similar to Anderson localization (fig. ITsFc)). A further increase in the impurity 
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(b) 
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Figure 18. Schematic localization scenarios (a) Pure bosonic superfluid in an 
optical lattice, (b) Shift of the effective potential depth due to fermionic impurities, 
(c) Localization by interfering paths of the bosonic wavefunction scattered by randomly 
distributed fermionic impurities, (d) Localization due to percolation. A random 
fermion distribution hampers the establishment of a coherent connection and causes 
the localization of bosonic ensembles in superfluid "islands". Copyright (2006) by the 
American Physical Society. 



density will lead to the formation of "forbidden walls". Once the walls in this quantum 
percolation scenario lead to a sufficiently complicated labyrinth like structure for the 
bosonic wavefunction, a single coherent superfluid phase can no longer be sustained and 



Definitely, further experiments are required to shine more light on the observed 
phenomena. A possible shift of the Mott-insulating transition can be probed by 
a spectroscopic investigation of the excitation spectrum of the mixture similar to 
probing of the excitation gap of the pure bosonic Mott insulator in [201 1128j . Lattice 
site occupation can be probed by rf spectroscopic techniques |129j . The correlation 
properties of the fermionic and bosonic components can be revealed by noise correlation 
techniques |130lll31ill32j . With the availability of heteronuclear Feshbach resonances, it 
is straightforward to extend our studies to almost arbitrary interaction strength between 
the fermionic and bosonic constituents, and scattering-length dependent measurements 
have recently been performed by Best et al. |133j . Also, the coherence properties of 
Bose-Bose mixtures with repulsive (interspecies-) interactions have been investigated by 
Catani et al. [134] . In this case, the discussion is complicated by the limited overlap of 
the clouds due to repulsive interactions; in fact, a significant shift of coherence reduction 
to lower lattice depths is observed even for close to vanishing overlap between the clouds. 
In any case, further studies will have to concentrate on isolation and identification of 
underlying phases and associated phase transitions in this rich system. 




Heteronuclear quantum gas mixtures 



43 



6. Ultracold heteronuclear molecules in a 3D optical lattice 

Being able to perform experiments with ultracold dense polar molecular gases is one of 
the outstanding challenges in atomic, molecular and optical physics |135j . The increasing 
interest in this problem is motivated by a wide range of of research directions which 
would directly benefit from the availability of such molecular samples. 

The quest for ultracold polar molecules 

A diatomic polar molecule is characterized by the presence of a permanent internal dipole 
moment occurring as a result of the electronic charge distribution not having the same 
"center of mass" as the nuclear charge distribution. The dipole moment associated with 
this electronic charge distribution gives rise to a long-range dipolar interaction between 
molecules. This interaction is novel in two senses: it is long range, and the interaction 
depends on the relative alignment of the constituents. In a sample of aligned dipoles, it 
is energetically more favorable for two dipoles to sit in front of each other than side by 
side. To date, all quantum gas experiments have studied effects of interactions which 
are short-range and isotropic, with the notable exception of experiments performed 
with ultracold Cr [136J which has a noticeable magnetic dipole moment. Studies of 
dipole-dipole interacting quantum gases are expected to reveal novel fascinating physics. 
For example, in the presence of dipolar interactions, the elementary excitations and 
the stability of the gas in various dimensions are expected to strongly depend on the 
dimensionality, the external confinement and the ratio of long range to short range 
interactions [HD [ll2]. To name only a few of the predicted effects, 

quantum hall states of dipolar gases have been considered, both for fermionic and bosonic 
particles |143l 1144] , phases of gases with dipolar interactions in optical lattices |140j and 
novel solitonic behavior [145] . 

It has been suggested by D. DeMille in 2001 [M] to use polar molecules for quantum 
computation. The idea is to encode quantum information in the orientation of the 
molecular dipole moment relative to an external field. Using the dipole alignment of a 
single molecule, a single qubit can be realized. A qubit register, i. e. multiple qubits 
can be realized by loading dipolar molecules into a ID optical lattice. In order to 
make this quantum information storage register a quantum computer, some degree of 
coupling needs to be introduced between the qubits. This is where the long range 
dipolar interaction becomes important: it couples molecules at different lattice sites and 
is used as the quantum bus. Following the initialization of the quantum register and 
the proper register operation, the readout of individual qubit states can be achieved by 
adding an electric field gradient along the direction of the optical lattice to separately 
address lattice sites through the Stark effect. Whereas trapped atomic ions are currently 
the most advanced system from an experimental point of view, the suggested quantum 
computation scheme is being evaluated as a promising candidate for future quantum 
computers and simulators. 

Paramagnetic polar molecules are currently being considered for precision 
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experiments testing upper limits on the permanent electric dipole moment of the 
electron (EDM). While the standard model sets an upper limit on such an EDM which 
is vanishingly small (less than 10~^°ecm), supersymmetric theories predict values of 
the EDM which are getting within reach of current experimental efforts. Improving 
the upper limit of the EDM would thus eventually enable certain new theories to be 
ruled out, while finding a non-zero EDM would be an exciting signature of T- violating 
physics beyond the standard model. It has initially been pointed out by Sandars |146] 
that the electronic EDM signature (a linear Stark effect) could be enhanced in heavy 
paramagnetic atoms, which is expressed by the so-called "enhancement factor". Several 
experiments are under way testing limits on the EDM using e. g. atomic fountains, 
atoms confined in solid helium matrices and beam experiments. In almost all atomic 
experiments, the achievable sensitivity is affected by the electric fields that can be 
applied without initiating a discharge. In paramagnetic heavy polar molecules, an 
enormous degree of polarization can be achieved due to the internal electric field with 
only moderate external fields. Experiments on polar molecules and the EDM are 
currently under way using e. g. PbO, TIF, and also molecular ions. Measurements 
of the EDM with polar molecules would greatly benefit from the availability of ultracold 
samples with long interrogation times, and the counting rate would be substantially 
increased by having all molecules in a well defined quantum state. 

Actively pursued routes to molecular degeneracy fall into two different categories, 
depending on whether the starting point is a deeply bound molecular sample (direct, or 
'molecular' approach) or whether molecules are assembled from precooled and possibly 
degenerate atoms (atomic or 'indirect' approach). 

The key obstacle in obtaining ultracold or even degenerate molecular samples 
compared to atomic Bose-Einstein Condensates or degenerate Fermi gases is the lack of 
laser cooling techniques for molecular gases. The main issue is the wealth of molecular 
energy levels, enhanced by rotational and vibrational degrees of freedom and the 
corresponding lack of clear cycling transitions for laser cooling of molecules. Currently 
followed routes to molecular degeneracy starting with deeply bound molecules (direct 
approach) therefore either rely on cryogenic techniques bringing the sample into contact 
with a thermal reservoir or on interaction with external electric fields slowing down the 
external degrees of freedom or filtering out low velocity classes. All of these can be 
followed by subsequent evaporative cooling of a trapped sample, but the temperature 
gap into the ultracold and even degenerate regime has yet to be bridged. Among the 
techniques being actively pursued are buffer gas cooling |147j . Stark deceleration [148] 
and velocity filtering |149j . 

The technical difficulties associated with slowing down the external degrees of 
freedom of molecules have made another route to molecular degeneracy attractive. 
This approach is based on ultracold atomic samples which are routinely produced in 
laboratories today. The goal is to assemble ultracold molecules in both their external 
and internal (rovibrational) ground state from ultracold atoms by means of a coherent 
photoassociation step |150[ 1151] . 
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The bottleneck in these photoassociation schemes will be the small Franck-Condon 
overlap between free atoms and deeply bound molecular states |15H 1152] . At this 
points, heteronuclear Feshbach molecule formation as demonstrated for the first time 
in this work |3] is a key step in obtaining a suitable initial state for coherent Raman 
photoassociation schemes. The Feshbach association provides us with a well-defined 
initial molecular state with known molecular distance and a much better Franck-Condon 
overlap with deeply bound molecular states. 

Recent years have seen a lot of progress in molecule formation based on Feshbach 
resonances. Molecule creation at Feshbach resonances is based on the fact that the 
occurrence of the resonance itself is due to the coupling between a bound molecular 
state energetically close to the incoming energy of two free atoms. 

When the energy of the bound molecular state is above the free-atom threshold, 
the coupling between the two close to degenerate states gives rise to strong attractive 
interactions between atoms. In free space, no two-body bound state exists near threshold 
(see fig. |8] right). 

When the energy of the bound molecular state is below the open channel, the 
interaction between atoms is large and repulsive, and pairs of atoms can occupy the 
two-body bound state and form a weakly bound molecule with typical binding energies 
between a few 10 kHz and a few MHz. Starting from an atomic sample, such molecules 
can e. g. be created by sweeping a magnetic field from attractive interactions through 
the resonance center position over to repulsive interactions (see fig. Islleft). 




In a groundbreaking experiment performed at JILA [84j, this magnetic field sweep 
technique was first used to create a coherent atom-molecule superposition. In this and 
following experiments with bosonic atoms and Feshbach resonances, it was observed 
that the resulting molecular lifetime was rather low and on the order of about 1 ms, 
limited by coUisional losses. A big step forward in Feshbach chemistry was achieved 
when it was realized that bosonic molecules created from two fermionic atoms using the 
magnetic field sweep technique would be coUisionally stable due to the Pauli exclusion 
principle |153j . The experimental demonstration of long-lived bosonic molecules created 
from fermionic atoms was soon to follow [881 EHl 11541 189| 192] . 

Parallel to this development, it was realized that the lifetime limitation of molecules 
created from bosonic atoms [8^ll55][86l[85| l87] can be overcome in 3-dimensional lattices 
by creating molecules from atom pairs in isolated wells of the lattice and thus inhibiting 
inelastic collisions [25l 1156] . Creation of molecules in optical lattices highlighted the role 
of the external trapping potential in molecule formation. A model system of two particles 
in a harmonic trap was already considered in 1997 |157] . and in 2005 experiments 
with optical lattices demonstrated that the existence of the external trapping potential 
shifts the free-atom threshold, resulting in the existence of two-body bound states on 
the attractive side of the Feshbach resonance where no stable molecules exist in free 
space jl58] . also seen later in the behavior of repulsively interacting pairs |159] . 

Today, many people consider molecule creation through magnetic field sweeps as 
the way to create molecules at Feshbach resonances. In recent years, however, a number 
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of techniques has been used to create Feshbach molecules. For example, it seems to be a 
unique property of ®Li that molecule formation can be achieved by evaporative cooling 
of thermal lithium atoms on the repulsive side of the Feshbach resonance without any 
magnetic field sweeps ^1] . Another technique which has been demonstrated is magnetic 
field modulation at a Feshbach resonance |160j . When the modulation frequency 
corresponds to the binding energy of the bound molecular state at the time average 
value of the magnetic field, free atoms can be converted into molecules on the repulsive 
side of a Feshbach resonance. Another possibility of creating molecules is rf association 
of atoms into molecules as developed within this work. 



Ultracold heteronuclear Feshbach molecules 

While heteronuclear Feshbach resonances had already been predicted in 2003 |161j (for 
— ^''Rb) and experimentally identified through coUisional losses in the middle of 2004 
[9^ [93] , all of the fascinating physics which has become accessible through Feshbach 
resonances has been limited to homonuclear systems. As a crucial step towards the 
exploration of dipole-dipole interacting molecular systems, novel quantum computation 
schemes and fundamental measurements, this work |4j presents the first experimental 
demonstration of ultracold long-lived heteronuclear molecules. Molecules are created 
within the single well of an optical lattice at a heteronuclear Feshbach resonance by 
means of a novel technique developed for this purpose and based on rf spectroscopy |162j . 
Lifetime and molecule creation efficiency have been measured and are consistent with 
a physical picture of lattice occupation. Both molecules stable in the absence of any 
external potential, attractively interacting atom pairs and repulsively interacting pairs 
with a positive "binding energy" have been identified. A detailed understanding of 
this model system has been developed based on a pseudopotential approach for the 
atomic interaction in cooperation with F. Deuretzbacher, K. Plassmeier, F. Werner, D. 
Pfannkuche and F. Werner |5]. The model consistently treats both the anharmonicity of 
the lattice potential and the general case of particles with unequal trapping frequencies; 
excellent agreement with numerical simulations has been obtained. 

As discussed above, various techniques have been employed for molecule production 
at Feshbach resonances in the literature. Within this work, rf association in an optical 
lattice as a novel technique for molecule production has been used. The idea of this 
association scheme is to start with a non Feshbach-resonant mixture in one spin state. 
In another spin state, a Feshbach resonance occurs at the same magnetic field, and 
the difference between this state's energy and the non-interacting limit in the same 
spin state is precisely the binding energy. By shining in an rf photon which provides 
the undisturbed energy corresponding to the spin fiip plus the additional interaction 
or binding energy, we can populate the two-body bound molecular state (see inset in 



fig. 20). In the presence of an optical lattice, molecule production through rf association 
provides us with a precise measurement of the binding energy at the same time, thus 
extracting a maximum of information about the energy spectrum. Rf association allows 
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Figure 19. "'"K 4^5*1/2 energy levels in the presence of an external magnetic field. 
The magnetically trapped 1 1 , | ) state and the ||,— 5) / states used for rf 

spectroscopy are highlighted in blue. 

molecule production without rapid magnetic field sweeps and uncertainties associated 
with magnetic field settling, thereby removing technical complications which can make 
a measurement of the binding energy quite cumbersome. 

In order to understand rf spectroscopy as used for molecule creation in this work, 
let us first look at the energy levels of ^^Rb and ^''K in the presence of a homogeneous 
external magnetic field. For J = I ±1/2, the well-known Breit-Rabi formula yields the 
energy levels in the presence of an external magnetic field. Figure [19] shows energy levels 
of ^'^K for fields of experimental interest. Due to the low hyperfine splitting, relatively low 
fields are required to drive "^^K into the Paschen-Back regime (compared to, for example, 
^''Rb). The resulting magnetic field insensitivity of transitions between neighboring 
levels at experimental fields is one of the reasons for the success of rf spectroscopy 
in this system. States are labelled both by their low-field {F = I + J) limit F and 
mp quantum numbers and by the independent m/ and mj quantum numbers in the 
Paschen-Back regime. 

The Feshbach resonance that was used for molecule production in this work occurs 
between ^''Rb atoms in the |1,1) and atoms in the |9/2, — 9/2) absolute ground 
states at a magnetic field of 547 G. The idea for determining the energy spectrum of the 
|9/2, —9/2) (8> |1, 1) state is to drive a transition between the "^^K |9/2, —7/2) state and 
the Feshbach-resonant |9/2, —9/2) state. The former features a scattering length which 
is independent of B over the magnetic field range studied in the experiment and small 
(given approximately by the background scattering length)]^ the energy of the latter 
strongly varies as a function of B across the resonance compared to the undisturbed 
|9/2,— 9/2) energy level. The interaction energy (binding energy) of the Feshbach- 
resonant state |9/2, —9/2) is thus given as the difference of the observed |9/2, —7/2) 

* The closest known heteronuclear Feshbach resonances in this state are located at 522 G and 584 G [61] 
and have theoretical widths of below 1 G. 
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|9/2, — 9/2) transition frequency from the undisturbed transition frequency which can 
be calculated from the Breit-Rabi formula or extracted from a measurement where no 



bosons causing the collisional shift are present (see inset in fig. 20). Note that this 
measurement of the collisional shift is up to the small constant offset caused by the 
non-resonant |9/2, —7/2) |1, 1) interactions. 

6.1. Experimental protocol 

The experimental implementation of this scheme is as follows: A quantum degenerate 
mixture of and ^^Rb in the "magic" optical dipole trap is prepared in the non 
Feshbach-resonant |9/2, —7/2) ® |1, 1) state using the exact same procedure as in our 
experiments on tuning of interactions. The magnetic field is set to a value close to the 
Feshbach resonance occurring at 547 G. In total, we allow 150 ms for the field to stabilize 
close to its final value. During the last 100 ms, the 3D optical lattice intensity is increased 
from zero to its final value using a linear ramp. At magnetic field values of 547 G, the 
'^"K |9/2, -7/2) |9/2,-9/2) transition occurs at about 80 MHz. In order to probe 
the energy spectrum for a given magnetic field, an rf pulse with a Gaussian amplitude 
envelope (1/e^ full linewidth of 400 /xs) and total duration of 800 /xs is irradiated using 
the standard evaporation antennajj] For a description of the versatile rf /microwave 



manipulation system used to generate these pulses, see Appendix A and |40l 1111] . 

About 1ms after the end of the pulse, the lattice intensity is ramped down using 
a linear ramp within 1 ms in order to reduce the kinetic energy, and the dipole trap is 
then suddenly switched off, marking the beginning of the time of flight. During time of 
flight, we detect both remaining ^'^K |9/2, —7/2) atoms and ^''Rb |1, 1) atoms as well as 
weakly bound molecules and ^''K atoms in the |9/2, —9/2) channel. 

At fixed magnetic field, the rf pulse frequency is varied in the vicinity of the 



undisturbed |9/2, —7/2) — > |9/2, —9/2) rf transition to obtain rf spectra as in fig. 20 
for 5=547.13 G. The spectrum shows the number of ^''K atoms in the |9/2, —9/2) state 
and weakly bound Feshbach molecules as a function of rf detuning from the pure 
|9/2,— 7/2) |9/2, — 9/2) transition at the set value of the magnetic field. The large 
peak at zero detuning occurs whether bosons are present in the optical lattice or not. 
It stems from lattice sites where only one fermion and no boson is present and is used 
as a precise magnetic field calibration in connection with the Breit-Rabi formulajftj 

The second peak occurring about 14 kHz above the undisturbed (atomic) peak 
is due to the existence of a Feshbach-resonant shift at lattice sites where one boson 
and one fermion are present. The strong interaction between the two atoms at the 
Feshbach resonance introduces a strong differential energy shift between the |9/2, —7/2) 
and the |9/2, — 9/2) levels. In this case, additional energy is required to drive the rf 

tt The pulse shape is chosen such as to avoid the side lobes in spectroscopy that usually result from 
probing a two-level system with square pulses. A detailed discussion of the observed pulse shapes can 
be found in liO] , 

ft Note that the observed 1/e^ half linewidth of 1.7 kHz is in excellent agreement with what we expect 
based on our Gaussian rf pulse [40^ . 
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Figure 20. Rf spectroscopy of "OK-^^Rb at 547.13 G in a A0.5E^^ deep optical 
lattice [4. Copyright (2006) by the American Physical Society. 



transition. This implies that in the presence of ^^Rb, the energy of the target state 
|9/2,— 9/2) is lower than in the reference level, which shows that the target state has 
a negative binding energy and that heteronuclear molecules have been formed through 
rf association. Forming heteronuclear molecules through rf association thus measures 
the binding energy of these molecules at the same time (about 14kHz in this case). In 



the particular case of fig. |20] recorded at a magnetic field of 547.13 G, the interaction 
between ^^K and '^^Rb is attractive which, in free space, would not allow molecules 
to be created because the bound molecular state is above threshold. As we shall see, 



the negative binding energy of the molecules in fig. 20 is due to the presence of the 
optical lattice which admits the presence of a two-body bound state even for attractive 
interactions. The resulting pair close to threshold will decay into free particles when the 
lattice potential is removed. In order to avoid potential confusion, we will use the term 
'attractively interacting atoms' for this type of atom pair in the following discussion. 

6.2. Energy spectrum 

In order to extract information about the energy spectrum over the whole magnetic field 



range of the Feshbach resonance, spectra as in fig. |20]have been recorded for a number of 
magnetic fields around 547 G. For each of the spectra, the molecular binding energy has 
been extracted by fitting both atomic and molecular peaks and taking the difference in 
position between the atomic and the molecular peak. The resulting spectrum is shown 



in fig. 21, recorded for two different values of the lattice depth, 40.5 E^^ and 27.5 E"^*^. 



The vertical bar in fig. |2T] represents the center position Bq of the Feshbach 
resonance as determined from ^°K-^^Rb interactions in an optical dipole trap (see our 
discussion on tuning of interactions). In the magnetic field range considered here, the 
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Figure 21. Observed energy spectrum of ''^°K *^Rb heteronuclear molecules and 
atom pairs at an individual lattice site for two different values of the lattice depth j4j . 
Copyright (2006) by the American Physical Society. 



interaction is attractive above Bq and repulsive below Bq. 

The energy spectrum shows two basic branches. The lower branch or molecular 
branch is characterized by the presence of a negative binding energy, and the upper 
branch, which is only observed for repulsive interactions, is characterized by a positive 
"binding energy". In a system of pure bosons confined in an optical lattice, the 
behavior of these repulsively interacting pairs has been studied in recent experiments in 
Innsbruck |159] . 

In the molecular branch, one can distinguish between attractive and repulsive 
interactions. Above Bq, we produce attractively interacting atom pairs |158l which 
exhibit a negative binding energy due to the presence of the lattice potential. The 
sample rf spectrum in fig. 20 has been recorded in this regime. As the resonance center 
is crossed, there is a smooth transition into "real" molecules which are stable even in 
the absence of an external potential. As the scattering length decreases from infinity 
(at the resonance) towards zero, these molecules become more and more deeply bound. 
The maximum binding energy observed in the experiment is about -130 kHz. 

Thus, using rf association, we have created heteronuclear molecules for the first 
time. The rf association process allows us to precisely determine the binding energy of 
the molecules while we create them. A key advantage of this scheme is that it provides 
a lot of information about the system without any potentially cumbersome magnetic 
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field ramps. 

Several aspects of this experiment have been crucial in obtaining these results. One 
is that when the rf is detuned from any resonance, the signal recorded in the target state 
|9/2, —9/2) is effectively zero. So the molecule signal which we measure in that channel 
is a signal on zero, and not a small dip in some fluctuating atom number. This was 
particularly helpful in the first stages of the experiment, where the molecule fraction 
was very small due to a small number of lattice sites occupied by one ^^Rb and one 
atom. 

A second aspect is imaging of the atoms and molecules. While in principle we 
could have used the Stern-Gerlach technique for state-selective detection in ''"K, this 
technique limits the time of flight to a certain minimum value given by the necessity 
to apply the field gradient between release from the trap and the actual imaging flash. 
This minimum time of flight can make it difficult to detect small numbers of molecules, 
in particular when in the optimization phase of the experiment, since after a large time 
of flight, the gas becomes very dilute. More importantly, with Stern-Gerlach imaging, 
the molecules cannot be imaged directly and need to be dissociated prior to imaging. 

Instead, we have chosen high field imaging in order to state-selectively image both 
'^''K atoms and molecules. High field imaging relies on the fact that at the high magnetic 
field close to the Feshbach resonance where the '^°K atoms are mostly in the Paschen- 
Back regime (Hyperfine structure of the Zeeman effect), it is possible to construct a 
cycling transition starting from any arbitrary magnetic sublevel of the ground state. In 
this case, the total shell angular momentum J, the nuclear spin I and the corresponding 
z projections mj and mj are good quantum numbers, and we have cycling transitions 
{rrij = ±1/2; m/ = — / . . . + J) — (mj = ±3/2; mj = —I . .. + /). For the case of 
weakly bound Feshbach molecules, one photon of the detection pulse will dissociate the 
molecule, and then subsequent scattering processes detect atoms in the corresponding 
spin state. This detection process is spin dependent because in general, the various 
optical cycling transitions for the different spin states differ by several 10 MHz at least. 
The particular scheme used in our experiments is described in Appendix D 



Crucial to rf spectroscopy and rf association as discussed in this work is the influence 
of environmental magnetic field noise and magnetic field drifts in combination with the 
magnetic field sensitivity of the rf transition used in the experiment. At the magnetic 
fields of interest for this work (547 G), the rf transition in "^^K used here has a small 
magnetic field sensitivity of 67kHz/G. This is mainly a result of ^^K being close to 
the Paschen-Back regime. The other key to the excellent quality of the rf spectra is 
suppression of magnetic field noise and drifts. By eliminating sources of AC magnetic 
noise, line triggering the experiment and using a custom built precision regulator, we 
have achieved a magnetic field reproducibility of 2.7 mG at 547 G (57 measurements on 
11 consecutive days, as seen from the atomic peak in rf spectroscopy), corresponding to 
a magnetic field reproducibility of ~ 5 ■ 10~^. 
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6. 3. Properties of heteronuclear Feshbach molecules in 3D optical lattices 

After discussing creation of heteronuclear molecules, let us look at some basic properties 
of these weakly bound dimers, such as a quantitative understanding of the energy 
spectrum, and the lifetime and creation efficiency of the molecules. These aspects are 
particularly relevant in view of possible applications, such as optical transfer into more 
deeply bound molecular states and ultimately the creation of a polar molecular quantum 
gas. The quantitative analysis presented here has been the result of a collaboration 
with F. Deuretzbacher, K. Plassmeier and D. Pfannkuche at Hamburg (1. Institut flir 
Theoretische Physik) and F. Werner (ENS, Paris) [5]. 

The simplest possible model for our system is that of two atoms, one and 
one ®^Rb, interacting via s-wave scattering where the scattering length as a function of 



magnetic field is given by equation (16), at an individual lattice site with the shape of 
the lattice potential approximated by a harmonic oscillator. This basic model has been 
brought forward by Busch and coworkers in 1998 |157j . In their article, the authors 
specifically consider interaction in the form of a regularized 5-potential between the two 
atoms and analyze the special case of two identical atoms with equal masses and trap 
frequencies. The corresponding Hamiltonian is shown to have an exact solution. 

In our case, there are several important differences from this simple model. The 
most obvious one is that we deal with a system where the two particles are not 
identical. The second is that for experimentally relevant lattice depths, the harmonic 
approximation for the lattice potential is only a very rough approximation — the 
anharmoni cities are significant, as we shall see. Another aspect, which we will not 
be concerned with here, are deviations from the regularized 6 potential in combination 



with equation (16). These deviations will be relevant for larger binding energies (see e. 
g. |103] ). but not for the regime of binding energies considered here. 

Our model Hamiltonian thus consists of the kinetic energy of the two atoms, of the 
trapping potential and the interaction potential: 

H = Ai A2 

2mi 2m2 

I 22,^ 22 

27rh'^a d ^ n 

H d{r)—r + Vanh{ri,r2) 

/i or 

where a is the s-wave scattering length between the two atoms, mi,2 the respective 
mass, fi = nil ■ m2/(mi + 777,2) the reduced mass, ri^2 the respective single particle 
coordinate, u!i^2 the harmonic trap frequency for atom 1 and 2 and Vanh contains 
anharmonic correction to the external trapping potential. Reference |157] proceeds 
by introducing relative and center of mass coordinates r = r 1 — r 2 and R = 
where M = mi + m2 is the total mass. By introducing frequencies 
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Figure 22. (a) Energy spectrum of a heteronuclear system (^'^K-^^Rb) for the 
pseudopotential model at a lattice depth of 40.5 E^"" [5 . (b) Experimentally measured 
binding energy (dots) compared to (dashed line) theory with literature parameters for 
width w and center position Bq of the Feshbach resonance and (red solid line) fitted 
curve resulting in = 546.669(24)sys((2)stat G [5]. Copyright (2008) by the American 
Physical Society. 



^rel ■■=\l 

Au := ■\JujI + , 



the resulting Hamiltonian takes the form 

H = - ^^CM + Imu'cmR' (29) 

= HcM + Hrel + Hco 

with one part Hcm involving only center of mass coordinates, one term Hrei containing 
only relative coordinates and a final term Hco potential coupling the two coordinates. 
The crucial point here is that in the absence of anharmonicities and for equal trap 
frequencies, the coupling term vanishes and the center of mass and relative motion 
completely decouple, allowing the identification of exact solutions. This is the special 
case treated by ref. |157j . 

This decoupling is absent in our general case of unidentical atoms trapped in 
an anharmonic lattice well and prevents us from directly using these exact solutions. 
Instead, given ref. |157j . we can write down exact solutions for the first two groups of 
terms Hcm and Hrei- Using a numerical procedure, the coupling term Hco is then exactly 
diagonalized in the basis set given by the eigenfunctions of the first two groups of terms, 
and we can obtain both wave functions and energy structure of the full Hamiltonian. 

Figure [22|^a) illustrates the results for the experimentally relevant case of a lattice 
depth of 40.5 -E^*^. The blue solid line shows the solution to the full Hamiltonian 



Heteronuclear quantum gas mixtures 



54 



as obtained by the exact diagonalization procedure, while the dashed hne shows the 
energy spectrum resulting from the non-coupling terms of ref. [157] alone. The difference 
between the two curves is most pronounced in the repulsively interacting atom branch of 
the spectrum, which is a result of the relative wave function extending over a significant 
part of the external potential, and the coupling between center of mass and relative 
motion is important. As the molecule size becomes smaller and smaller for larger 
negative binding energies, the coupling term Hco becomes less important. 

Figure [22|^b) compares the solution to the full coupled problem to the experimentally 
measured binding energies. As can be seen, there is a good agreement between our 
experimentally measured binding energy and theory; however, there seems to be a 
slight shift between experiment and theory which is most likely due to insufficient prior 
knowledge in the determination of the parameters of eq. ( 16 ) connecting the scattering 
length and the magnetic field. By fitting the theoretical result for the binding energy to 
the experimental data (red solid line), we can thus obtain a new value for the Feshbach 
resonance center position Bq = 546.669(24:) systi'2') stat G. Details of the fitting procedure 
and the error analysis can be found in ref. [5]. 

In view of possible applications, the lifetime of weakly bound Feshbach dimers 
has been a particular concern since their first demonstration |8l]. In general, these 
molecules are very fragile objects due to their highly excited internal state, leading to 
large coUisional losses, possibly both molecule-molecule and atom-molecule collisions. 
Initial experiments with molecules created from bosonic atoms |8H I155i |86| |85| |87] 
correspondingly showed a very short lifetime, meaning that degeneracy could be 
achieved, but not thermal equilibrium. This limitation has been overcome in the 
presence of a deep 3D optical lattice [251 1156j . Molecules created from two fermionic 
atoms have shown a long lifetime close to Feshbach resonances due to Pauli-forbidden 
inelastic decay |153j . For heteronuclear molecules composed of a bosonic and a fermionic 
atom, the situation is somewhat more complicated; suppression of collisions due to their 
fermionic character is expected for more deeply bound molecules [93j- In addition, 
collisions with the remaining free atoms are expected to strongly limit the lifetime of the 
molecules, and this effect will tend to become more important the deeper the molecule 
is bound. 

In the experiment described here, heteronuclear molecules are produced within 3D 
optical lattices where the large tunneling time strongly reduces inelastic collisional losses. 
Figure [23|b) shows a measurement of the lifetime of the heteronuclear molecular sample 
in an optical lattice with a depth of 40.5 Ef'^ as a function of magnetic field. The 
lifetimes are the 1/e times obtained from exponential fits to the molecule number decay 
curves. We find a lifetime of about 120 ms for weakly bound attractively interacting 
atom pairs. In the vicinity of the resonance, the lifetime is about 80 ms, and drops to 
20 to 40 ms for more deeply bound molecules. 

So far, we have excluded tunneling between neighboring lattice sites from our 
discussion. Tunneling is not relevant for the short timescales involved in the molecule 
creation process. The presence of the optical lattice results in a very long tunneling time 
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Figure 23. (a) Cartoon picture of expected lattice occupation [5 . (b) Lifetime of the 
molecular sample as a function of magnetic field [5]. Copyright (2008) by the American 
Physical Society. 



between 10 and 20 ms for ^°K. The tunneling time for ^^Rb is larger as a result of the 
larger mass. Ultimately, however, tunneling becomes important, as we shall see. For a 
cartoon picture of our expected lattice occupation, see fig. [23|^a). In the optical lattice, 
prior to molecule creation, we expect a global occupation of at most one atom in 
the F = 9/2, rnp = —7/2 state per lattice site (Pauli exclusion principle). As far as 
^''Rb is concerned, we expect a central occupation number of three to five bosons per 
lattice site, dropping to one in the outer regions of the optical lattice. It is in this outer 
shell, at lattice sites with one '"^K and one *^Rb atom, that the rf association creates 
molecules. Once the molecular association has happened, the sites with molecules are 
surrounded by other sites with atoms that are quantum-mechanically distinguishable 
from the molecules (^°K in the F = 9/2,mF = —7/2 state and ®^Rb). There is some rate 
for tunneling of these remaining atoms to the lattice sites occupied by the molecules. 
Once the tunneling process has happened, the molecule and the atom can undergo 
three-body recombination, and the probability for this process can be determined from 
our simple theoretical model. Following Petrov et al. |153] , the decay probabihty (after 
tunneling has happened) can be calculated as the probability of finding all three particles 
(the two constituents of the molecule, plus the atom that has tunneled to the molecular 
site) in a very small volume. For the tunneling atom, we use the wave function of a non- 
interacting atom trapped at a lattice site. For the molecule, we use the wave function 
obtained from our solution to the coupled Hamiltonian (29 ). The resulting theory curve. 



scaled by a global factor to allow comparison to the experiment, is shown in fig. |23[a) 
and is in good agreement with the experimentally observed magnetic field dependence 
of the lifetime. Ultimately, this means that while the molecular lifetime in the optical 
lattice is already excellent for applications (such as transfer to lower lying molecular 
states and ultimately the absolute ground state), the lifetime can be further increased 
by removing any unpaired atoms from the optical lattice. At high magnetic fields, this 
can be achieved by a combination of rf pulses and resonant light beams. After completion 
of this work, one of us (S. O.) has joined a project aimed at transferring heteronuclear 



Heteronuclear quantum gas mixtures 



56 




attractively interacting 
atoms 



1.2 

1 

0.8 



0.2 



(b) 




expt. data 






theory — 




f t 






repulsively interacting 


atoms 





(a) 5-iC-3 54C.4 .54C.5 54G.C 54C.7 546.8 54C.9 547 

B I G 



(b) 



545. C 545.8 



54C 546.2 

B I Q 



546.4 546.6 



Figure 24. Rf association efficiency for botli tlic molecular (a) and the repulsively 
interacting atom pair brancli (b). Copyright (2008) by the American Physical Society. 



molecules created in dipole traps into the absolute internal molecular ground state. This 
project has published an extensive study of stability of molecules in the dipole trap 
when they collide with various atomic collision partners |163j . 

An important aspect of molecule production as demonstrated in this experiment 
is the efficiency of rf association, in particular as a function of magnetic field. There 
are two different aspects to the overall efficiency. A fundamental limit on the achievable 
absolute efficiency in any current molecule creation experiment is imposed by the number 
of doubly occupied lattice sites as compared to the overall atom number. This relation 
can be somewhat optimized experimentally by playing with temperature and overall 
particle numbers. 

A second aspect is the relative efficiency as a function of magnetic field and rf 
pulse power (i. e. at which fraction of the lattice sites occupied by one ^^Rb and one 

atom does the rf pulse create molecules?). In our experiments, in order to find 
suitable parameters for the rf pulses, we have driven Rabi oscillations on the atomic 
transition and then set the parameters such as to always obtain a tt pulse on the atomic 
transition. For all of our experiments, we have preserved this choice of parameters. 
Figure 24 shows the resulting relative transfer efficiency as a function of magnetic field. 
We expect that for heteronuclear scattering lengths in the final state which are close 
to the initial state {F = 9/2, mp = —7/2), i. e. at magnetic fields significantly above 
the resonance center or otherwise far detuned, we should create attractively interacting 
atom pairs at all available lattice sites occupied by one *^Rb and one ^'^K atom. The 
experimental transfer efficiency as determined from the ratio between the atomic and 
molecular peak in the rf spectroscopy has therefore been normalized to one in this region. 

In |5], we have established a model for the association process based on a Rabi 
model. After the rotating wave approximation and after integrating out spatial degrees 
of freedom, we obtain Rabi flopping on the molecular association transition with a 
Rabi frequency which is reduced by a factor of {^i\^f) compared to the atomic Rabi 
frequency. Here, is the overlap integral between the initial and flnal motional 
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states of the rf association process as determined from the pseudopotential model. 

Since we have chosen constant pulse parameters in the experiment, we expect the 
transfer efficiency to drop as the initial and final states become more and more dissimilar. 
This can be seen very clearly in fig. 24, both in the experimental data and in the 
prediction of the Rabi model. The quantitative agreement, however, is not complete. 
A further possibility of testing the hypothesis of the Rabi model is to look for Rabi 
oscillations between atoms and molecules on the molecular rf transition. This could 
provide further insight into the association process. 



7. Conclusions and Outlook 

During our two theses [39l HOj, we have realized an experimental apparatus for the 
preparation of quantum degenerate Fermi-Bose mixtures in 3D optical lattices with 
tunable interactions. We have presented important milestones towards a mixed statistics 
many-body system as a versatile model system for the simulation of quantum many- 
body Hamiltonians and demonstrated important steps on the way towards quantum 
degenerate polar molecular gases with long-range dipolar interactions. 

Starting from the realization of the so far largest particle numbers in a magnetically 
trapped degenerate ^°K-^^Rb system, we have investigated interaction effects due to 
the large and attractive background interaction between the components. Tuning of 
heteronuclear interactions in the vicinity of a Feshbach resonance has been demonstrated 
for the first time within this work and has allowed studies of arbitrary interactions 
between bosons and fermions. The complete phase diagram of harmonically trapped 
mixtures has been accessed. Starting with experiments on stable density distributions 
for attractively and repulsively interacting mixtures, for strong repulsive interactions, 
phase separation has been observed and a Feshbach-induced collapse for strong attractive 
interactions. A heteronuclear p-wave resonance has been identified as described in [2], 
confirming a theoretical assignment of Feshbach resonances in the ^°K-*^Rb system and 
opening intriguing perspectives for tuning of the anisotropy of the interaction. These 
studies on heteronuclear Feshbach resonances pave the way for studies with strongly 
interacting Fermi-Bose mixtures. 

From the very beginning, the experimental setup described here has been designed 
for studies in 3D optical lattices. Within this work, Fermi-Bose mixtures have been 
loaded into a three-dimensional optical lattice for the first time, demonstrating a large 
fermion-induced loss of coherence of the bosonic sample. The results are currently the 
subject of intense theoretical analysis, and explanations ranging from disorder-related 
localization phenomena over adiabatic thermodynamic effects and a mean-field induced 
shift of the superfluid to Mott-insulator transition are being discussed. Further insight 
into these observations may be gained by implementing additional diagnostics. This may 
include Bragg spectroscopy, excitation spectroscopy to probe possible excitation gaps in 
deep optical lattices, probing of lattice occupation through rf or microwave spectroscopy 
and noise correlation analysis. The studies may also be extended to mixtures with 
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varying heteronuclear interaction strength (see most recent results by Best et ai, |133j ). 

In ID or 2D optical lattices, one or several spatial degrees of freedom can 
experimentally be frozen out, resulting in the preparation of low-dimensional systems. 
Depending on interaction strength and lattice occupation, charge-density wave or pairing 
phases may be observable. Low dimensional systems are of particular interest for direct 
comparison of experimental results to theory. For 1-dimensional situations (2D optical 
lattice) exactly solvable models exist and direct comparison to quasi-exact DMRG 
calculations is possible. 

The combination of tunable interactions with 3-dimensional optical lattices has 
resulted in the first demonstration of ultracold heteronuclear Feshbach molecule 
formation. Long-lived heteronuclear Feshbach molecules have been created in a 
controlled fashion through rf association of atoms to molecules. This has allowed a 
precise determination of the binding energy of the molecular sample as a function of the 
magnetic field. Lifetimes between 20 ms and 120 ms have been experimentally observed, 
and excellent agreement has been found between experimentally determined properties 
of the molecules (binding energy, lifetime) and a a simple universal theoretical model. 

Molecule formation may be the starting point for studies of pairing phases in the 
optical lattice, including studies of 3-body bound states. After molecule association and 
removal of the left-over atomic fraction, the lattice potential may be ramped down. The 
molecule formed from a bosonic ^''Rb and a fermionic atom is again a fermion, and 
the observation of a molecular Fermi sea may be demonstrated. 

In the optical lattice, Feshbach molecules are created in their absolute external 
ground state at a lattice site, but in a highly excited rovibrational internal state. Pulsed 
two-color photoassociation may be used to transfer these molecules into their internal 
ground state. The resulting molecule would exhibit a permanent electric dipole moment 
(polar molecule). The dipole moment gives rise to a long-range, anisotropic interaction. 
In the optical lattice, this long-range interaction may be used to implement quantum 
computation schemes. In the many-body limit, these polar molecules give access to 
novel quantum gases with long-range anisotropic interactions. 

After completion of this work, one of us (S. O.) has joined a project |163l 1103] 
aimed at optically transferring jl64j weakly bound Feshbach molecules into the absolute 
ro-vibrational and electronic ground state, and most recently realized a dense ultracold 
all ground state polar molecular gas |165j . 

Appendix A. Rf / microwave manipulation 

In our experiments, we have used a novel rf / microwave manipulation scheme |11H BO]. 
where all timinig, frequency and amplitude control is derived from one fast arbitrary 
frequency generator. The versatility and simplicity of operation as well as the precise 
timing down to the nanosecond level has greatly enhanced the experimental flexibility. 
The setup was developed based on the observation that many of the commercially 
available frequency generators provide a plethora of modulation options, whereas they 
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Figure Al. rf/microwave manipulation setup. Amplitude, frequency and timing 
control is provided by the VFG-150. Its output is connected to an rf multiplexer (M), 
which either directly feeds the signal to an amplifier and a coil-shaped antenna (up 
to 150 MHz), or to two mixers ((8)), summing the output with two fixed frequencies. 
Undesired sidebands in the resulting signal are eliminated with filters (F), and the 
signal is amplified and sent through a circulator to eliminate back reflections from the 
antenna. The resulting power is delivered to the atoms using a coil shaped device (for 
'*'^K hyperfinc manipulation close to 1.2 GHz) or using an open end waveguide (for 
^^Rb hyperfinc manipulation at 6.8 GHz). Adapted with permission from llllf - 



are usually unable to rapidly and accurately perform for example simple linear frequency 
sweeps or AM modulation with a given envelope. The usual Adiabatic Rapid Passage 
(ARP) sequences for Zeeman state preparation require sweep times on the order of 10 ms, 
whereas the dwell time of many commercial synthesizers, i. e. the time duration of each 
individual "slot" in a linear sweep, is limited to precisely 10 ms; the usual workaround is 
to apply an analog voltage sweep to the FM modulation input. Also, many synthesizers 
need a lot of time for frequency switching both between frequencies and from fixed 
frequency operation to FM, sometimes as long as a second. Often, each operation 
requires its own synthesizer, resulting in an expensive and complex 'synthesizer stack'. 

In this experiment, we have implemented a solution which relies on just one 
synthesizer (VFG-150, available through Toptica Photonics, and developped by T. 
Hannemann in the group of C. Wunderlich), offering precise control of amplitude, 
frequency and phase every 5 ns. The synthesizer covers an output frequency range 
of several 100 kHz to 150 MHz. The device (equipped with a USB2.0 input) is remote 
controlled from a small server application running on a PC. We have made the source 
code of this application available at http : / /vfg- control . sourcef orge . net/ Many of 
the required rf manipulations can be performed directly using the output of the VFG- 
150 (evaporation, state preparation at low magnetic field and rf spectroscopy of ''^K 
at high magnetic fields). For other applications, such as transfer from one hyperfine 
manifold to the other, the required frequencies exceed the direct output range of the 
VFG-150. To cover all these frequencies, we have connected the output of the VFG-150 
to an rf multiplexer which either directs it directly to an amplifier and antenna or to 
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a series of mixers which mix it with fixed frequency microwave sources. We have used 
mixing with a 1.2 GHz fixed frequency to generate sweeps and pulses for ^"K hyperfine 
manipulation and with a 6.8 GHz fixed frequency to transfer '^''Rb from F = 2 to F = 1 
(see fig. Al). 

The advantage of this approach is that just one synthesizer is being used. The 
overall setup is versatile to operate; it allows all kinds of shaped pulses (e. g. the 
Gaussian amplitude envelope pulses used for molecule creation) and offers a 5 ns timing 
accuracy. At the same time, this scheme is highly cost-effective because it requires only 
one versatile source (the VFG-150) and inexpensive fixed frequency sources. 



Appendix B. 'Magic' optical dipole trap 

Compared to magnetic trapping, optical dipole traps offer the advantage of confining 
an arbitrary spin state. Furthermore, the value of the magnetic field becomes an 
experimental degree of freedeom. This capability is important for experiments with 
Feshbach resonaces as described in this article. A special interest is in relatively weak 
confinement, allowing optical lattice geometries which are as close as possible to a 
homogeneous situation. 

However, when working with shallow traps and mixtures, the differential 
gravitational sag of the two components becomes an important issue. An atom with 
mass m confined in a harmonic trap with trap frequency u in the presence of gravity 
experiences a gravitational sag given by zq = —{g/u'^). This expression is found by 
calculating the local minimum of the combined gravitational and trapping potential. 
Different atomic species and even the different spin states of the same atomic species will 
in general experience a different trapping frequency when confined in the same magnetic 
or optical dipole trap, giving rise to a difference in gravitational sag. The importance 
of this issue depends on the value of the trap frequency. It is generally a non-issue 
for the compressed magnetic trap or for strong dipole traps. It does however play a 
significant role in the weak isotropic magnetic trap used for transfer of laser precooled 
atoms and for experiments in shallow and weak dipole traps such as the ones used in 
experiments with optical lattices aiming at a situation close to homogeneous systems. 
The difference in gravitational sag can then lead to a significantly reduced overlap of the 
two species. For example, for a far detuned optical dipole trap (equal trapping potential 
for both species) with a vertical trap frequency of 50 Hz, the differential gravitational 
sag is 49 /im! 

The "magic" crossed dipole trap discussed below has allowed us to successfully 
overcome this limitation: if we can make the trap frequency in the direction of gravity 
the same for both species, the gravitational sag will automatically be the same for 
^''Rb and "^"K. In an optical dipole trap, this can be done by a suitable choice of the 
dipole laser wavelength for confinement in the direction of gravity. This so-called "magic 
wavelength" can be calculated to be 806. 7 nm taking into account contributions from 
the Dl and D2 lines. 
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Figure Bl. Measurement of the magic wavelength for our dipole trap. 



With this idea in mind, we first performed a measurement of the magic wavelength 
by measuring the ratio of and ^^Rb trap frequencies in a dipole trap as a function 



of wavelength and found the magic wavelength to be 807.9 nm (see fig. Bl). Based on 
this information, a crossed beam optical dipole trap was designed with one beam at 
the magic wavelength providing confinement in the z (vertical) and x direction. The 
potential minimum of the atoms in this beam is made to overlap with the position in 
the magnetic trap. We tend to hold the intensity of this beam constant. A second beam 
is derived from our Yb:YAG lattice laser (A = 1030 nm) and provides confinement in 
the y and z direction. This beam is overlapped with the atomic position in the "magic" 
beam. This is also the beam that we use for optical evaporation in the crossed dipole 
trap by lowering its intensity after transfer of atoms into the optical trap. This optical 
dipole trap allows us to go to weak optical dipole traps {u = 271 ■ 50 Hz) without having 
issues with reduced overlap of the two clouds. It has been the basis for our experiments 
on tuning of interactions, mixtures in lattices and heteronuclear molecules. A detailed 
discussion of this novel dipole trap can be found in [39] . 



Appendix C. Optical lattice 

Part of our experiment is a 3D optical lattice setup. The setup uses an Yb:YAG disc 
laser with an output power of 20 W as a laser source. In order to avoid motional heating 
in the optical lattice due to a 'shaking' interference pattern, both frequency and intensity 
of the laser have been stabilized. The short term frequency stabilization locks the laser 
frequency to a reference cavity and has required the custom development of frequency 
control elements for the laser cavity. On longer time scales, the reference cavity follows 
the laser frequency. The achieved short term linewidth is 20 kHz relative to the reference 
cavity. Details can be found in |110[ HQ] . 



Figure CI shows a general view of the beam preparation scheme for the optical 
lattice setup. Behind an optical diode and pickups for laser power monitoring and 
frequency stabilization, the light from the lattice laser is split up into 4 independent 
beam paths using polarization optics. Each of these beams passes through an Acousto- 
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Figure CI. Optical lattice setup with beam preparation, intensity stabilization, 
coUimation unit and retroreflection. Not shown: Imaging along lattice axes and edge 
filters for overlapping with MOT beams. 



Optical Modulator (AOM) in single-pass for intensity control and is coupled into an 
optical fiber. All of the AOM frequencies are detuned by roughly 10 MHz with respect 
to each other in order to make the interference terms between the individual beams 
so fast that they are averaged out on any experimentally relevant time scale. Three 
of the beam paths are used for the three orthogonal retrorefiected lattice beams with 
mutually orthogonal polarizations, and one of them is used for one of the two beams of 
the "magic" crossed optical dipole trap. 

After passing through the optical fiber, each of the lattice beams is coUimated and 
focused. A pickup is used for the intensity stabilization. The focussing lens has a 
focal length of 400 mm. Having passed this lens, the lattice beam is overlapped with 
the 3D MOT beams using dichroic mirrors and delivered to the atoms. The standing 
wave configuration for the optical lattice is produced by separating the MOT and 
lattice beams again after passing the glass cell using an identical "disentangling" dichroic 
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mirror, recoUimating the beam using a second lens identical to the focussing lens and 
retroreflecting it. Not shown in fig. CI are the dipole trap beams and imaging systems 
along all of the three lattice axes. These allow us to overlap the lattice beam focus with 
the atoms and to optimize the lattice alignment. This is done by nulling out oscillations 
that result from a quick switch-on of the lattice in the presence of lateral misalignment. 



Appendix D. Imaging 

In our experiment, we detect atoms using resonant absorption imaging. The standard 
absorption imaging procedure is to record three images, one where the optical density 
of the atoms attenuates the beam falling onto the CCD (absorption image), one where 
the laser beam alone is imaged (reference image) and one image without any light 
(dark image). From these three images, the optical column density of the cloud can 
be calculated (see e. g. |166j ). We have used several variations of this scheme in our 
experiments. All of these rely on properties of modern CCD chips allowing the user to 
record two images in very rapid succession (interline transfer CCD chips, such as the 
SONY ICX205AL). Typically, the interframing time is on the order of /is, if not less, 
and after two images, the readout takes on the order of hundreds of ms, larger than a 
useful time of fiight in absorption imaging. We have made use of this capability to take 
two images in very rapid succession during time of fiight in two different ways: 



Two-species imaging 

In a first scenario (two-species imaging), we have used the first double image to store 
absorption images of and ^^Rb, respectively. The second double image contains the 
reference image for both species, and the last double image contains the dark image for 
both. This allows us to record both species on one camera and on one imaging axis in 
the same experimental run. This capability has been crucial in the initial optimization 
of the experiment and in the exeriments on collapse dynamics. 



Fast imaging 

A second scenario makes use of the capability to take two images in rapid succession to 
eliminate the typical interference fringes from the images. Both the absorption and the 
reference image typically contain interference fringes caused by the coherent nature of 
the imaging light and etaloning at various surfaces in the optical path. These fringes 
are removed in the standard normalization procedure as long as they are stationary 
over experimental time scales. Usually, reference and absorption image are separated 
by the readout time of the CCD chip, thereby making the result sensitive to drifts of 
the interference pattern on acoustic time scales. Using the two images taken with /is 
separation as absorption and reference image, the interference patterns are stationary 
in between the two images, and we obtain a much better image quality and sensitivity 
to small atom numbers. Since the atoms do not fall very far in the short interframing 
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time, wc need to make sure the imaging beam does not 'see' them when the reference 
image is taken. For ^^Rb, which wc have typically detected out of the F = 1 state when 
this imaging procedure was used, we take the reference image first, when the imaging 
beam resonant with the F = 2 ^ F = 3 transition on the D2 line does not see the 
atoms. We then shine in a repumping pulse which transfers the atoms to F = 2 and 
take the absorption image after the reference image. 

For ^°K, we have taken advantage of the fact that we have usually detected at high 
magnetic fields (see below). This allows us to take the absorption image first at high 
magnetic fields, then quickly switch the 547 G magnetic field off in a couple of 10 //s, 
which drives the atoms out of resonance with the imaging beam and lets us take the 
reference image. The improved image quality and sensitivity to small atom numbers 
has been crucial to the measurements on tuning of interactions, mixtures in 3D lattices 
and heteronuclear molecules. 

State-selective imaging 

We have used two different techniques of state-selective imaging to identify the 
distribution of the atoms over the various sublcvcls of the hyperfine ground state 
manifold. The first and more common technique is Stern-Gerlach imaging, which 
uses a strong magnetic field gradient applied during time of fiight to spatially separate 
the spin states. Stern-Gerlach imaging has the advantage of simultaneously providing 
information over all spin states. However, it requires the application of the field gradient 
for several milliseconds, thereby introducing a minimum time of fiight which can impede 
the detection of small atom numbers. A different technique of state selective imaging, in 
particular when working at high magnetic fields, e. g. when using Feshbach resonances, is 
to take advantage of the transition from the Zeeman to the Paschen-Back regime. In this 
regime, the detection transitions mj = ±1/2 to mj = ±3/2, with Am/ = 0, are cycling 
transitions, and we can build a cycling transition from any sublevel in the ground state 
hyperfine manifold. Detection at high magnetic fields has the advantage of allowing 
a short time of fiight, and allowing the detection of both atoms and heteronuclear 
molecules. However, it requires the atoms to be sufficiently far in the Paschen-Back 
regime, and for each spin state to be imaged, a separate detection laser and CCD camera 
is required. This technique has been crucial to our measurements on heteronuclear 
molecule formation. 

In our particular system, the detection transition for the Feshbach-resonant state 
F = 9/2, mp = -9/2 (mj = -1/2, mj = -4 ^ nij = -3/2, mi = -4) and for 
weakly bound molecules is detuned by -764 MHz from the standard zero field cycling 
transition (F = 9/2, = 9/2 ^ F = 11/2, = 11/2). The closest other state, 
which is F = 9/2, = —7/2 at low field and rrij — — 1/2, mj — —3 at high field, has 
a cycling transition which is detuned by -833 MHz from the aforementioned standard 
zero field transition. Hence, we have a frequency difference of 69 MHz between the 
detection transitions for the two states used in the rf spectroscopy, allowing us to state- 
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selectively image both weakly interacting atoms in the F = 9/2, mp = —7/2 channel 
and Feshbach molecules and strongly interacting atom pairs in the F = 9/2, rrip = —9/2 
channel. Note that in our particular system, high field detection is only possible with 
the component because at the magnetic fields of interest, ^''Rb is closer to the low 
field (Zeeman) regime than to the Paschen-Back regime, and the corresponding cycling 
transitions are not available. 

In order to perform the detection, we leave the Feshbach magnetic field on during 
time of fiight and detection. For imaging, we use two beams derived from our cooling 
and detection laser and differing by 69 MHz in frequency. After laser cooling and during 
the evaporative cooling, we detune our imaging laser by about 700 MHz |167j to make 
it resonant with the high field detection transition. This high field detection technique 
has allowed us to state-selectively image both atoms and weakly bound molecules with 
very short time of fiight, high sensitivity and unparalleled simplicity (no dissociation 
of molecules or ramping of fields). Another beneficial aspect of this detection scheme 
is that in combination with fast imaging (see above), it has allowed us to drastically 
reduce the typical interference fringes in absorption imaging, again improving on the 
signal to noise. Details of this scheme can be found in [ID]. 
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